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ABSTRACT 


Static  constitutive  stress-strain  relations  are  developed  for  concrete 
at  intermediate  pressure  levels,  up  to  10-12  ksi  mean  normal  stress.  An 
elastic-plastic  model  is  developed  that  qualitatively  fits  the  loading 
features  exhibited  for  a  variety  cf  controlled  laboratory  tests  conducted 
in  this  program  and  found  in  the  literature.  Laboratory  tests  were 
conducted  on  several  batches  of  concrete  with  3/8-inch  maximum  aggregate. 
Considerable  attention  was  paid  to  the  casting  of  the  concrete  so  that 
batch-to-batch  repeatability  could  be  obtained,  and  to  specimen  preparation. 
Some  data  are  presented  for  aggregate  variation  from  3/16-inch  to  3/4-inch 
so  that  scaling  to  different  aggregate  could  be  performed.  The  laboratory 
tests  presented  represent  a  series  of  tests  where  all  stresses  and  strains 
have  been  measured,  thereby  allowing  the  shear  and  dilatation  stress-strain 
responses  to  be  observed  during  a  variety  of  load-unload  paths. 

Strain  rate  effect s  are  not  considered  in  the  model  developed,  nor 
are  the  features  of  unloading  correctly  handled.  Each  of  these  represents 
an  area  where  additional  detailed  work  is  needed. 
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Concrete  is  a  composite  material  of  naturally  occurring  sand  and  gravels 
or  crushed  rock,  bonded  together  by  some  form  of  cement,  and  may  contain  as 
many  as  eight  different  components.  These  Include  coarse  aggregate,  sand, 
unhydrated  cement,  cement  gel,  gel  pores,  capillary  pores,  air  voids  and 
free  water.  For  many  engineering  purposes, concrete  may  be  considered  as  a 
composite  of  aggregate,  cement  matrix,  water  and  unfilled  voids.  This 
simplified  model  does  not  In  any  way  imply  that  all  concretes  are  similar  In 
material  response.  Quite  the  contrary,  'concrete'  alone  is  not  aoequate  to 
specify  a  given  material.  The  nature  of  concrete,  including  the  complicated 
chemical  reactions  that  take  place,  is  well  documented,  for  example  by 
K.  Newman  (Ref.  1). 

The  response  of  concrete  to  either  very  rapid,  quasistatic  or  long-term 
loading  Is  quite  different,  and  sufficient  data  are  available  to  Indicate 
that  concrete  exhibits  a  relatively  complicated  stress-strain  response,  very 
different  from  the  frequently  used  approximation  of  linear  elasticity.  Under 
quasistatic  loading  nonlinear  effects  become  particularly  significant  at 
stresses  greater  than  about  one-half  the  maximum  stress,  as  indicated  for 
example  by  the  papers  shown  In  Reference  2,  "The  Structure  of  Concrete"  (1968). 
Under  very  rapid  loading,  concrete  shows  viscous  response,  and  a  recent  report 
by  Read  and  Malden  (Ref.  3)  gives  a  good  survey  of  the  response  of  concrete 
to  shock  loading. 

Even  though  plain  concrete  Is  a  composite  material,  the  work  Indicated 
here  treats  It  as  a  homogeneous  continuum.  That  is,  the  aggregate,  cement 
matrix,  water  and  voids  are  assumed  to  be  distributed  In  a  homogeneous  manner. 
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and  the  geometry  of  the  concrete  specfnen  is  assumed  sufficiently  large  so 
that  the  size  of  the  aggregate  and  voids  are  relatively  small  compared  to 
the  geometry  of  the  specimen.  Testing  methods  have  recently  become  available 
to  determine  the  bulk,  or  continuum,  response  of  concrete  for  a  variety  of 
loading  conditions  (Ref.  4),  which  thereby  makes  the  continuum  approach  to 
modeling  the  stress-strain  response  particularly  attractive.  Admittedly 
for  some  applications,  such  as  shock  loading  where  the  load  pulse  may  be 
small  compared  with  size  of  the  aggregate  and  voids,  the  continuum  model  may 
need  additional  detail.  A  reinforced  concrete  structure  would  need  to  be 
treated  as  a  'composite'  structure  consisting  of  the  reinforcement  in  a 
concrete  matrix,  and  the  continuum  model  developed  here  would  suffice  to 
model  the  matrix. 

The  work  here  considers  only  quasistatic  loading  and  does  not  address 
creep  or  high  strain  rates.  The  objective  of  this  research  was  to  develop 
general  three-dimensional  constitutive  equations  which  qualitatively  fit  the 
features  of  deformation  exhibited  by  concrete  at  intermediate  pressures,  up 
to  10-12  k s 1 ,  and  to  this  end  a  variety  of  loadings  were  made  on  specimens 
from  a  specific  laboratory  cast  batch  of  concrete. 

A  constitutive  model  was  developed  that  represented  the  data  obtained 
in  this  program.  Although  this  model  fits  only  the  specific  concrete  studied, 
it  could  be  applied  to  another  concrete  of  interest  by  using  different  values 
for  the  constants  involved. 
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Strength :  Considerable  attention  has  been  given  to  the  tensile  and  compressive 
strength  of  plain  concrete.  The  effect  of  mix  ratios,  aggregate  size  and 
grading,  curing  tine  and  environment,  have  generally  been  Investigated  (Ref.  1 
and  5).  The  water-cement  mix  ratio  appears  to  have  the  greatest  influence 
on  strength,  with  lower  water  ratios  giving  greater  strength.  Different 
mixing  and  casting  techniques  can  give  different  amounts  of  entrapped  air, 
and  hence,  large  differences  in  porosity. 

Specimen  design,  fabrication  and  end  conditions,  or  constraints,  have 
likewise  been  investigated.  Although  'casting'  of  specimens  to  the  desired 
size  and  shape  has  been  extensively  used  (Refs.  1,5,6),  the  surface  mortar 
layer  and  nonuniformity  of  specimens  due  to  the  different  casts  may  be  less 
desirable  than  casting  a  single  large  block  and  diamond  cutting  or  coring 
samples  from  the  single  casting.  Compressive  specimens  with  length  to 
diameter  ratio  of  two  are  generally  considered  adequate  to  obtain  a  region 
with  uniform  stress,  while  shorter  specimens  will  undoubtedly  exhibit  non- 
homogeneous  stress  through  the  specimen  (Ref.  7).  The  type  of  end  lubrication 
used,  if  any,  seems  to  vary;  however,  tests  on  hard  rocks  by  Wawersik  (Ref.  8) 
suggest  that  'lubricants'  may  lead  to  low  apparent  strength  caused  by  intrusion 
In  pores  and  splitting  of  the  end.  End  conditions  which  lead  to  apparent 
increased  or  decreased  strength  must  carefully  be  separated  from  the  effect 
of  triaxlal  loading  on  strength.  Although  end  effects  due  to  unknown  end 
conditions  lead  to  triaxlal  stress  loading  in  the  specimen,  the  results  will 
likely  not  be  scalable  or  subject  to  interpretation. 

Tests  under  triaxlal  compressive  stress  (Refs.  9-13)  show  that  the  compressive 
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strength  of  concrete  is  sensitive  to  confining  pressure  and  Increases 
rapidly  with  increases  in  confining  pressure.  Figure  1,  taken  from  Newnan 
and  Newman  (Ref.  5),  illustrates  this  effect.  This  Figure  also  shows  that 
the  classic  expression  of  Richart,  et  al.  (Ref.  9) 

oj  =  fc  +  4.1  o, 

is  inaccurate  at  the  higher  confining  pressures. 

Concrete  tensile  strength  is  much  lower  than  compressive  strength,  but 
the  ratio  of  these  strengths  depends  on  the  mix  parameters  (Ref.  14).  A 
ratio  of  uniaxial  compressive  strength  to  tensile  strength  of  about  ten  may 
be  considered  typical  although  ratios  on  the  order  of  twenty  have  been 
observed  (Ref.  13). 

Tests  under  triaxial  extensio  stress  (Refs.  12,13)  and  under  biaxial 
stress  (Refs.  15,16,17)  indicate  that  strength  is  not.  only  sensitive  to 
confining  pressure  but  is  also  more  precisely  dependent  upon  thp  "tress 
state.  For  example,  the  maximum  stress  difference  is  great'  triaxial 

compression  than  in  triaxial  extension  at  the  same  confining  pressure,  as 
shown  by  Figure  2  which  is  a  replot  of  data  by  Hobbs  (Ref.  13). 

The  tension  plus  biaxial  compression  data  of  McHenry  and  Kami  (Ref.  18), 
Eresler  and  Pister  (Ref.  19)  and  others  as  summarized  by  Newman  and  Newman 
(Ref.  5)  all  indicate  that  a  steady  reduction  in  tensile  strength  results  as 
the  biaxial  compressive  stress  is  increased.  This  can  be  seen  in  Figure  2 
which  also  indicates  that  as  the  lateral  compressive  stress  is  raised,  a  point 
is  reached  where  the  axial  stress  is  zero.  This  should  correspond  to  the  biaxial 
strength  under  equal  biaxial  stresses,  and  thus  provides  a  comparison  point 
between  the  two  kinds  of  tests.  Further  increases  in  the  lateral  compressive 
stress  produce  failure  under  axial  compressive  stress,  again  as  illustrated 
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Figure  1.  Effect  of  confining  pressure  on  strength 
of  concrete,  from  Newman  and  Newman  (5). 
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gure  2.  effect  of  triaxial  stress  state  on  strength  of  concrete 
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in  Figure  2. 

The  strength  of  concrete  under  biaxial  compressive  loading  has  been 
investigated  by  Kupfer,  Hilsdorf  and  Rusch  (Ref.  15)  using  special  brush 
platens  to  minimize  adverse  end  conditions.  A  strength  increase  of  16  percent 
in  equal  biaxial  compression  relative  to  uniaxial  compression  was  obtained, 
which  is  in  accord  with  the  value  of  10  to  25  percent  suggested  by  Newman 
and  Newman  (Ref.  5).  This  value  Is  not  universally  obtained;  for  example, 
the  extension  tests  of  Chinn  and  Zimmerman  (Ref.  12)  show  an  increase  of 
approximately  110  percent  in  equal  biaxial  compressive  strength  over 
uniaxial  compressive  strength. 

A  number  of  failure  criteria  have  been  proposed  to  predict  concrete 
strength  under  general  mul tiaxial  stress  conditions..  A  recent  review  by 
Newman  and  Newman  (Ref.  5)  indicates  that  this  is  not  a  settled  issue.  The 
classic  Mohr-Coulomb  criteria  states  that  strength  can  be  described  by 

°1  '  °3  =  f(°l  +  °3> 

This  expression  has  been  used  by  a  number  of  investigators.  For  example, 
the  expression  by  Richart,  et  al .  (Ref.  9)  given  previously  can  be  put  in 
this  form  with  f  being  a  linear  function.  However,  the  Coulomb-Mohr  hypothesis 
implies  that  the  strength  value  does  not  depend  on  the  intermediate  principal 
stress,  in  variance  with  biaxial  test  results.  Bresler  and  Pister  (Ref.  19) 
and  McHenry  and  Kami  (Ref.  18)  have  used  a  generalization  of  the  Mohr-Coulomb 
criterion  in  which  the  octahedral  shear  stress  is  taken  as  a  function  of  the 
octahedral  normal  stress.  A  good  fit  to  data  in  the  mixed  tension-compression 
stress  range  was  obtained. 

Chinn  and  Zimmerman  (Ref.  12)  foun  '  that  neither  the  Mohr  theory  nor  the 
octahedral  stress  generalization  fit  thjir  compression  and  extension  test 
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results.  This  conclusion  was  supported  by  Mills  and  Zimmerman  (Ref.  20). 

They  proposed  that  their  data  could  be  fitted  by  using  a  modification  of 
the  octahedral  shear  stress,  however,  the  application  of  this  theory  was 
not  presented. 

Stress-Strain  Response:  During  loading  in  uniaxial  compression ,  concrete 
behaves  in  a  nearly  linear  elastic  manner  up  to  about  half  the  maximum  stress, 
where  significant  microcracking  begins.  Unloading  from  above  this  stress 
level  will  result  in  permanent  set,  while  continued  loading  will  eventually 
result  in  massive  break  up  of  the  concrete  microstructure  and  large  scale 
dilation  as  maximum  stress  is  neared  (Refs.  1,5).  Careful  control  of  the 
testing  machine  strain  rate  allows  the  complete  stress-strain  curve  to  be 
obtained  (Refs.  8,21),  and  shows  that  fracture  does  not  occur  immediately 
after  reaching  peak  load.  The  load  tends  to  drop  with  increasing  strain, 
with  large  scale  cracking,  slabbing  and  complete  loss  of  cohesion  occurring 
well  after  peak  load.  Specimens  loaded  to  slightly  beyond  peak  load  and 
then  unloaded  may  show  no  macrocracks,  and  in  general  have  the  appearance 
of  an  uiideformed  specimen. 

The  stress  level  at  the  onset  of  significant  mi crocra eking  is  termed 
the  "discontinuit1'  stress"  (Ref.  5)  and  can  be  observed  by  deviations  from 
linearity  in  tie  stress- strain  response.  The  microcracking  has  been  detected 
by  a  number  Oi  methods  including  microscopic  and  X-ray  examination  of  specimen 
exterior  and  cross  section  surfaces,  measurement  of  acoustic  emissions,  and 
other  methods  (Kef.  5).  A  measurement  of  all  three  principal  strains  permits 
the  calculation  of  the  relative  change  in  volume.  Increases  in  volume  (relative 
to  the  elastic  deformation)  due  to  microcracking  have  been  observed  by  a 
number  of  investigators  (Ref.  5). 

Concrete  is  capable  of  large  deformations  under  confining  pressure. 
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Chinn  and  Zimmerman  (Ref.  12}  applied  axial  strain  of  16%  to  triaxial  com¬ 
pression  specimens  at  confining  pressures  of  25  ksi  or  more,  and  many  other 
investigators  have  noted  similar  increases  in  maximum  axial  deformation 
(Refs.  10,11,13).  The  stress-strain  response  becomes  significantly  non¬ 
linear  at  hig:.  deformations  and  permanent  set  upon  unloading  is  large. 

Measurement  of  all  principal  strains  has  not  received  widespread 
ttention.and, in  general,  insufficient  data  are  available  to  define  the 
triaxial  stress-strain  response.  Also  relatively  little  da  Lu  I  s  t  to 
determine  the  unloading  and  reloading  path,  the  effect  of  varying  amounts 
of  free  water  (dry  or  saturated),  temperature  and  thermal  cycling  or  other 
preconditioning  effects  on  the  triaxial  stress-strain  response. 

Concrete  Constitutive  laws:  Little  previous  work  has  been  performed  on 
developing  constitutive  equations  that  even  reasonably  represent  the  triaxial 
stress-strain  response  of  concrete.  This  is  likely  due  to  two  factors; 
first,  that  experimental  techniques  for  measuring  all  the  independent  com¬ 
ponents  of  the  strain  tensor  with  sufficient  accuracy  have  only  recentlv 
become  available,  ana  secondly,  the  usage  of  large  computers  that  can 
effectively  utilize  more  accurate  (and  perhaps  complicated)  stress-strain 
representations  has  only  recently  become  common.  Because  of  this,  the  stress- 
strain  behavior  of  plain  concrete  has  usually  been  modeled  as  either  elastic 
or  elastic-perfectly  plastic.  Nilson  (Ref.  22)  has  presented  a  typical  computer 
program  for  structural  analysis  of  reinforced  concrete  using  nonlinear 
properties.  However,  an  adequate  representation  of  the  necessary  properties 
was  not  available. 

On  the  basis  of  a  careful  study  of  the  stress-strain  behavior  of  several 
rocks,  Swanson  (Ref.  23)  and  Brown  and  Swanson  (Ref.  4)  proposed  that  a  strain 
hardening  plasticity  model  could  be  used  to  phenomenologically  represent  the 
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features  of  dilatancy  (bulking)  and  permanent  set.  Baron,  et  al.  (fief.  24) 
have  added  to  the  model  a  strain  hardening  "capped"  yield  surface  that  was 
particularly  appropriate  for  modeling  soils,  because  the  capped  yield 
surface  permitted  some  control  over  the  amount  of  dilatancy  exhibited  by 
the  plasticity  model.  Recently  Swanson  (Ref.  25)  has  shown  how  his  earlier 
strain  hardening  model  could  be  combined  with  a  capped  yield  surface  in  a 
nur.-.er  appropriate  for  some  porous  rocks. 

Herrmann  (Ref.  26)  has  recently  given  a  survey  of  constitutive  equations 
for  the  shock  compaction  of  porous  materials,  but  only  briefly  touches  on 
general  three-dimensional,  stress-strain  response. 
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SECTION  III 


EXPERIMENTAL  TECHNIQUES 

Concrete  Mix:  Smooth  aggregate  was  selected  which  contained  very  small 
percentages  of  crushed,  irregular,  or  elongated  material.  The  primary 
mineralogical  components  are  sandstone,  quartzite,  limestone/marble, 
granite,  gneiss,  and  schist,  which  are  naturally  occurring  around  Salt  Lake 
City,  Utah.  The  aggregate  is  generally  rounded  to  subangular  because  it 
is  stream  deposited  material.  To  achieve  a  minimum  aggregate  to  specimen 
cross  section  ratio  of  one  to  five,  3/8  inch  was  selected  as  the  maximum 
aggregate  size  for  the  detailed  tests  on  thick-walled, hollow  cylinders  and 
solid  cylinders,  with  some  unconfined  compression  tests  run  ori  different 
batches  of  concrete  having  3/16,  1/2,  and  3/4  inch  maximum  aggregate  sizes 
to  investigate  the  strength  dept  -ence  on  maximum  aggregate  size. 

Grading  curves  were  selected  for  the  four  different  aggregates  so  as 
to  maintain  the  overall  surface  area  per  unit  of  concrete  weight  (in2/lb) 
constant  within  the  resolution  of  the  calculations,  estimated  at  +30%. 

This  was  considered  important  since  surface  area  is  a  major  factor  in 
establishing  the  water/ cement  ratio.  The  surface  area  was  calculated  for 
the  aggregate  by  assuming  all  aggregate  was  spherical  in  shape.  This 
calculation  was  thus  only  approximate. 

The  aggregate/cement  ratio  by  weight  of  6.3/1  was  selected  to  afford 
a  'medium  workability'  (Ref.  27)  of  the  3/8-inch  aggregate  mix.  Figure  3 
shows  the  actual  grading  curves,  where  the  percentages  are  given  as  weight 
percentages  of  the  total  aggregate  passing  through  the  screen  size  indicated. 

A  single  lot  of  aggregate  was  obtained  and  stored  for  all  batches  cast,  and 
a  laboratory  check  analysis  was  run  on  each  aggregate  mixture  to  determine  the 
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Figure  3.  Actual  grading  curves  for  batches  cast 


actual  grading  curve. 

A  water/cement  ratio  (by  weight)  of  0.6,  and  a  Type  I  Portland  cement 
were  selected.  To  ensure  homogeneity  of  the  cement,  a  local  supplier  ran 
a  special  batch  which  was  bagged  directly  from  the  mixer.  This  one  lot  was 
used  for  all  concrete  batches  cast.  Table  I  presents  the  batches  cast. 

Table  I 

Concrete  Specimen  Batches 


Batch  No, 

Maximum 

Aggregate  Size, 

1 nches 

Comments 

1 

3/3 

1  tub,  inadequate  vibration 

2 

3/8 

1  tub 

3 

3/4 

1  tub,  different  surface  are e 

4 

3/16 

1  tub 

5 

3/4 

1  tub 

6 

1/2 

1  tub 

7 

3/8 

1  tub,  plus  5  hollow  cylinders 

B 

3/8 

5  hollow  cyl inders 

Mixing,  Pouring  and  Curing:  The  desired  amount  of  graded  aggregate  was 
weighed  and  placed  in  a  5-cubic  foot  drum  mixer,  and  mixed  several 
minutes  to  ensure  an  even  distribution.  A  sample  was  taken  for  moisture 
analysis,  and  the  mixer  was  sealed  with  a  clear  plastic  film.  A  moisture 
analysis,  which  entailed  careful  weighing,  then  heating  to  250"F  and  weighing 
at  30-minute  intervals  until  nc  further  weight  loss  was  observed,  was  run 
on  the  aggregate  sample.  A  dry  aggregate  weight  of  1%  was  allowed  as  normal 
interstitial  water  content,  and  the  remainder  was  considered  as  excess  water 
which  reduced  amount  of  water  added  to  t he  mix. 

Appropriate  cement  and  water  quantities  were  weighed  and  added  to  the 
mixer  to  bring  the  water,  cement,  IX  moisture  aggregate  ratio  by  weight  to 
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0.6  /  1  /  6.3.  The  batch  was  mixed  5  minutes  and  then  poured  into  approximately 
2-foot  diameter  tubs  (to  be  cored  later  for  solid  cylinders)  or  in  forms 
for  hollow  cylinders.  The  tubs  were  vibrated  15-20  seconds  with  a  1%-inch 
diameter  commercial  concrete  vibrator  and  the  forms  for  thick  walled  cylinders 
lesser  amounts.  In  general,  this  amount  of  vibration  was  adequate  to  stop 
air  bubbles  from  appearing  on  the  surface  and  to  produce  only  very  slight 
amounts  of  water  at  the  surface. 

The  concrete  was  moved  immediately  to  a  constant  temperature-humidity 
curing  room.  After  4  days  the  forms  were  removed  and  the  cast  specimens 
covered  with  saturated  burlap  material  kept  moist  by  absorption  f;om  water 
reservoirs.  Solid  cylinder  specimens  were  wet  diamond  core  drilled  from 
the  tubs  at  about  20  days  and  then  placed  under  water  for  continued  curing 
(except  for  a  few  specimens  as  iioted  later).  Some  specimens  were  removed 
from  the  water  at  24  days  to  be  ground  on  sides  and  ends,  strain  gaged,  and 
then  tested  at  28  days.  Other  specimens  remained  submerged  until  about 
80  days  and  then  were  kept  in  the  constant  temperature  (78°F)/constant 
humidity  (o0%)  room  until  time  of  test. 

Specimen  Preparation  -  Solid  Cylinders:  Solid  cylinders  were  core  drilled  to 
about  2  3/4-inch  diameter  using  a  water  lubricated  diamond  core  bit  on  a 
large  vertical  drill  press  as  shown  in  Figure  4  with  each  specimen  marked  for 
identification  and  location  in  the  larger  block.  By  using  proper  core  bit, 
feed,  speed,  lubrication,  and  very  rigid  alignment,  solid  cylinders  were 
obtained  with  no  appreciable  chipping  or  loosening  of  the  exposed  aggregate. 

No  indication  of  internal  cracking  was  suggested  by  either  microscopic 
observation  or  by  sc.  iter  In  tensile  or  compressive  strengths.  The  solid 
cylinders  were  ground  to  2.7-inch  diameter  with  a  tolerance  of  .001  inch 
on  roundness  and  straightness.  Thereafter  they  were  cut  and  ground  to  length, 
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(6  inches  for  most  of  the  tests)  parallel  and  flat  to  within  .0003  inch. 

Specimens  for  unconfined  tests  had  strain  gages  applied  directly  to  the 
concrete  surface  with  commercial  strain  gage  adhesives.  Typically  gages 
were  applied  at  the  specimen  center  in  axial  and  transverse  orientations 
spaced  180°  apart  around  the  specimen.  Foil  gages  with  a  3/4~inch  gage 
length  (Micro  Measurements  EA-06-750DT-120)  were  selected  for  most  of  the 
tests,  as  discussed  later.  Terminal  strips  and  leads  were  attached,  and  the 
specimens  appeared  as  shown  in  Figure  5a.  For  tensile  testing,  aluminum 
alloy  tension  grips  were  epoxy  bonded  to  the  specimen  ends,  as  shown  in 
Figure  5b. 

Specimens  for  tests  under  confining  pressure  had  large  surface  voids 
filled  with  a  fine  commercial  grout,  and  then  a  .005 -inch -thick  soft 
aluminum  sleeve  was  shrunk  onto  the  specimen  with  about  a  .005-inch  inter¬ 
ference  fit.  Strain  gages  were  bonded  to  the  aluminum  jacket  as  shown  in 
Figure  rc.  and  the  specimen  was  then  sealed  in  tygon  tubing  as  shown  in 
Figure  5d.  Tn-*  effect  of  the  aluminum  sleeve  on  strength  was  considered 
negligible  (Refs.  4,28). 

S pen' me  Preparation  -  fallow  Cylinders:  The  thick-wall ed,  hoi  1  ow -cyl i nder 
speciiiie.s  t)  be  tested  in  axial  compression  and  internal  pressurization  (hoop 
tension)  were  cast  in  forms  7  inches  inside  diameter  by  11  inches  outside 
diameter  by  20*  inches  long.  The  metal  forms,  as  illustrated  in  Figure  6a, 
were  removed  af+er  4  days, and  the  specimens  were  covered  inside  and  out  with 
saturated  buriap  material.  At  about  20  days  the  specimens  were  placed  under 
water  until  about  80  days.  The  specimen  ends  were  diamond  ground  flat  and 
perpendicular  to  the  specimen  axis,  to  within  .001  inch.  The  inside  of  the 
cylinders  were  lightly  sandblasted  with  fine  grit  to  expose  subsurface  voids 
(which  might  collapse  under  internal  pressure),  and  these  voids  were  filled 
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C.  CONFINED  COMPRESSION 
OR  TENSION 


d.  COMPLETELY  SEALED  SPECIMEN 


Figure  5 


Compression  and  tension  test  specimens. 


18 


with  a  fine  commercial  grout-  Because  of  the  matrix  cement  surface  layer 
inherent  in  cast  specimens,  2-inch-long  foil  strain  gages  (Micro  Measurements 
EA-06-20CBW-120)  were  selected  (as  opposed  to  3/4 -inch -long  gages}  for  use 
on  the  hollow  cylinders.  The  gages  were  applied  at  about  the  specimen 
center  internally  and  externally  at  axial  and  transverse  directions  180° 
apart.  A  specimen  ready  for  testing  is  shown  in  Figure  6b. 

Testing  -  Solid  Specimens:  Specimens  were  tested  using  a  very  stiff,  well 
aligned  servocontrolled  press.  The  press  for  unconfined  tension  and  com¬ 
pression  tests  is  shown  in  Figure  7.  Load  was  measured  from  a  precalibrated 
load  cell  and  strains  from  high  elongation  strain  gages  directly  mounted 
on  the  test  specimen  as  previously  discussed.  Axial  strain  rate  was  kept 
constant  at  10' V second  by  the  servocor.trol  1  ed  press,  and  because  of  the 
high  stiffness  of  the  press  the  compressive  stress-strain  curve  beyond 
maximum  stress  was  recorded  in  some  cases. 

For  tests  on  solid  cylinders  under  confining  pressure,  a  test  vessel 
and  loading  press  shown  in  Figure  8a  were  used.  The  specimen  was  aligned 
inside  the  test  vessel,  with  a  precal ibrated  load  cell  and  pressure  measuring 
manganin  coil  inside  the  test  vessel.  Strains  were  again  measured  by  high 
elongation  strain  gages  directly  mounted  on  the  specimen,  with  a  pressure 
correction  of  +0.5  x  10'8/p$i  (Ref.  29)  applied  to  all  strain  gage  readings 
(which  was  generally  insignificant  for  the  pressures  tested  here).  Kerosene 
was  used  as  the  fluid  inside  the  test  vessel,  with  pressure  generated  by  a 
servocontrolled  intensifier  shown  in  Figure  8b.  The  manganin  pressure  coil 
inside  the  test  vessel,  a  strain  gage  on  the  specimen,  or  the  load  cell, 
provided  the  input  signal  to  control  the  intensifier  unit  such  that  the  desired 
pressure  or  load,  or  the  desired  specimen  strain  condition  was  met.  Axial 
strain  rate  was  generally  kept  at  10"Vsecond,  as  for  the  unconfined  tests. 
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All  signals  were  recorded  in  analog  form  with  direct  signal  calibration. 

Load  and  pressure  are  accurate  to  ±2%,  while  strains  are  somewhat  less  accurate 
with  uncertainty  mainly  due  to  the  porous,  inhomogeneous  nature  of  the  concrete. 
Repeatability  of  stresses  and  strains  from  similar  tests  on  adjacent  concrete 
specimens  was  ±3-4%.  No  specimen  end  lubrication  was  used  for  compression 
tests;  the  concrete  was  loaded  directly  against  a  steel  end  platen. 

Preliminary  Tests:  The  sign  convention  adopted  throughout  this  report  is 
that  tensile  stress  and  extensional  strain  are  positive,  and  compressive 
stress  and  shortening  strain  are  negative.  Axial  refers  to  length  of  the 
specimen  and  transverse  to  the  diameter.  No  apparent  anisotropy  was  noted 
from  visual  or  microscopic  examinations,  and  no  difference  in  strength  or 
stiffness  was  noted  for  different  directions;  therefore,  no  distinction  has 
been  made  between  transverse  strains,  and  in  general,  isotropic  behavior 
is  assumed. 

Variation  in  Strain  Gage  Size:  Previous  work  reported  in  the  literature  has 
suggested  that  strain  gage  readings  on  concrete  and  mortar  specimens  are  a 
measure  of  tne  average  strains  to  within  at  least  5%,  as  long  as  the  gage 
length  •  greater  than  four  times  the  maximum  aggregate  size  (Ref.  30).  To 
better  rsurminr  the  effect  of  strain  gage  size  on  strain  reading,  several 
2.7-inch  diameter  by  6.0-inch  long  specimens  were  prepared  and  tested  in 
unconfined  compression  with  2.0-inch  long  (>5  times  maximum  aggregate  dimension), 
0.75-inch  long  and  0. 5-inch  long  strain  gages  applied  directly  to  the  specimen 
and  adjacent  to  each  other,  to  indicate  axial  strain.  The  strain  indication 
difference  between  the  longest  and  shortest  gages  was  less  than  8%,  with 
essentially  no  difference  between  the  2.0-inch  and  0.75-inch  long  gages,  as 
shown  in  Figure  9a.  Since  experience  had  revealed  that  larger  strain  gages 
are  more  subject  to  failure  at  confining  pressure  because  they  are  more  likely 
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to  cover  and  be  deformed  into  specimen  voids,  the  gage  length  selected  was 
0.75  inch,  which  is  twice  the  maximum  aggregate  size  (for  most  tests). 

Cylinder  Length/Diameter  Ratio;  Newman  and  l.achance  (Ref.  31)  and  Johnson 
and  Sidwell  (Ref.  7)  suggest  a  length  over  diameter  (1/d)  ratio  of  > 2.5  to 
ensure  freedom  from  end  effects  in  the  central  gage  section.  To  verify  an 
adequate  length  to  diameter  ratio  for  our  test  techniques,  several  2.7-inch 
diameter  by  6.0-inch  long  specimens  were  instrumented  with  three  0.75-inch 
long  strain  gages  to  indicate  transverse  strains  over  about  a  2-inch  center 
section  length  of  the  specimen.  These  were  specimens  from  batch  2  cured  about 
28  days.  One  gage  was  located  at  the  center  of  the  specimen  and  one  gage 
1  inch  above  and  1  inch  below  center,  with  lateral  strain  selected  as  the  most 
sensitive  indication  of  end  effects  since  typical  specimen  barreling,  due 
to  the  end  constraint,  would  be  indicated  by  lateral  gages.  The  three 
lateral  strains  were  uniform  to  <5/  difference  for  specimens  tested  without 
lubrication  between  platens  and  specimen,  as  shown  in  Figure  9b. 

To  further  investigate  the  effect  of  specimen  configuration  on  apparent 
physical  properties,  unconfined  compression  tests  were  run  on  specimens  with 
1/d  ratios  of  1.1,  2.2  (similar  to  the  previous  tests)  and  3.3  (nominal 
3,  6,  and  9-inch  lengths  by  2.7-inch  diameter).  These  were  specimens  from 
batch  2  cured  for  about  90  days,  but  not  under  water.  The  stress-strain 
curves  indicate  that  all  specimens  exhibit  an  initial  tanget  modulus  (slope 
of  axial  stress-axial  strain  curve  at  near  zero  stress)  of  4.7xl06  psi  ±8* 
as  seen  in  Figure  10a;  however,  maximum  loads  were  971C,  827 D,  and  7660  psi 
for  1/d  ratios  of  1.1,  2.2,  and  3.3,  respectively,  as  plotted  in  Figure  10b. 
This  represents  a  17/  apparent  strength  reduction  from  3  to  6-inch  specimen 
length,  but  only  a  7%  reduction  from  G  to  9  inches,  further  indicating  that 
a  1/d  of  2.2  is  adequate  for  tests  here  {7%  is  considered  small,  of  the 
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Figure  10.  Compressive  loading  of  specimens  with  different 
length-to-diameter  ratios. 


order  of  batch  scatter  as  will  be  shown  later).  From  these  tests  it  was 


(  .  v  *nd  effects  are  negligible  for  our  techniques,  if  a  1/d  ratio 

d. 

Van  n  Aggregate  Size:  To  determine  the  effect  of  maximum  aggregate 

size  on  strength,  28 -day  unconfined  compressive  tests  were  performed  on 
specimens  with  maximum  aggregate  sizes  of  3/16,  3/8,  1/2,  and  3/4  inch 
(2.7-inch  diameter  by  6.0-inch  long  specimens).  Aggregate  specific  surface 
area  was  held  constant  by  selecting  grading  curves  to  give  an  appropriate 
distribution  of  various  sizes  of  aggregate  in  each  mix.  Other  concrete 
parameters  were  also  held  constant,  as  listed  below. 

Cement 

Water-cement  ratio 
Aggregate-cement  ratio 
Mixing  time  and  technique 
Pouring  time  and  technique 
Vibrating  time  and  technique 
Curing  time  and  technique 
Specimen  preparation 

Maximum  stress  increases  only  slightly  with  increasing  aggregate  size  as 
shown  in  Figure  11.  A  strength  decrease  of  less  than  113  was  observed 
with  decrease  in  maximum  aggregate  size  from  3/4  to  3/16  inch.  These  values 
cannot  be  compared  directly  with  other  published  tests,  for  example  that  given 
by  Walker  and  Bloem  (Ref,  32)  for  3/8  and  3/4-inch  aggregates,  since  no  attempt 
wus  made  to  hold  specific  surface  area  constant  for  different  aggregate 
nixes  in  these  other  testing  programs. 

One  concrete  hatch  with  3/4-inch  maximum  aggregate  (batch  3)  was  mixed 
with  a  specific  surface  area  (inVlb)  about  20-50%  lower  than  the  other  batches. 
The  strength  was  17%  lower  than  the  other  3/4-lnch  aggregate,  batch  5,  with  the 
weakening  believed  to  occur  because  the  lower  surface  area  tics  up  less  water 
in  surface  bending  and  c'fcctively  increases  the  water  available  for  hydration, 
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Figure  11.  28-day  compressive  strength  of  various  batches  cast. 
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or  equivalently  this  gives  an  effective  increase  In  the  water-cement  ratio, 
which  is  known  to  produce  weaker  concrete  (Ref.  27).  In  another  batch,  about 
50%  less  vibration  during  casting  resulted  in  excessive  air  entrapment  and 
substantially  weakened  concrete,  as  shown  by  the  poStrt  for  baUh  1  on 
Figure  11. 

From  these  tests  on  variable  aggregate  size  concretes  it  is  concluded 
that  the  strength  Is  more  dependent  on  aggregate  surface  area  than  on 
aggregate  size,  and  that  careful  aggregate  grading  with  consideration  for 
specific  surface  area  is  required  to  produce  concrete  with  similar  properties. 
Furthermore,  such  sensitivity  to  pouring  and  casting  processes  points  out 
the  need  for  extreme  quality  control  in  obtaining  test  specimens. 
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SECTION  IV 


EXPERIMENTAL  OATA 

All  tests  reported  in  this  section  were  conducted  on  2.7-inch  diameter 
by  6.0-inch  long  solid  cylinders,  cast  and  aged  for  greater  than  200  days. 

The  tests  were  performed  on  servocontrol led  machines  that  maintained  a 
constant  axial  strain  rate  of  10*4/second.  Specimens  were  nominally  'dry', 
in  that  they  had  been  exposed  to  dry  laboratory  environment  for  several  days 
prior  to  test,  and  tests  were  performed  at  room  temperature,  about  70°F. 

Average  density  at  the  time  of  test  was  2.27  +  .03  g/cc. 

Unconfi ned  Tests :  Figure  12a  shows  the  results  of  compression  tests,  with 
engineering  stress  plotted  versus  engineering  strain,  while  Figure  12b 
shows  the  tension  test  results.  Each  curve  is  the  average  of  2-4  tests  with 
material  scatter  shown  as  ba-s  on  the  curves.  The  apparent  elastic  modulus 
is  about  6x10°  psi  in  compression  and  7xl0c  psi  in  tension,  with  an  unconfined 
compressive  strength  of  7020  psi  and  a  tensile  strength  of  805  psi. 

If  the  test  machine  is  adequately  control  Table,  the  decreasing  stress 
part  of  the  stress-strain  curve  can  be  readily  mapped  as  was  done  in  Figure  13 
where  a  compression  specimen  was  loaded  beyond  maximum  stress  and  unloaded, 
then  reloaded  until  large  scale  slabing  occurred  (as  shown  by  the  picture  in 
insert).  Visible,  prominent  axial  cracks  and  slabing  occur  well  after  maximum 
compressive  stress.  In  tension  fracture  occurs  mainly  through  the  mortar, 
with  fracture  seldom  running  through  the  aggregate  (Refs.  1,5).  This  fracture 
pattern  Is  quite  different  than  for  extension  tests,  which  are  discussed  later. 

The  aging  of  concrete  is  generally  well  documented  (Refs.  1,5)  and  shows  a 
rapid  increase  in  strength  during  the  first  28  days  followed  by  slower  increase 
in  strength  for  the  next  60  days,  and  relatively  little  increase  in  strength 
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after  90-100  days.  The  strength  of  the  concrete  tested  here  is  plotted 
versus  log  time  in  Figure  14. 

Constant  Confining  Pressure  Tests:  Compression  tests  were  conducted  at 
constant  confining  pressure  of  0.5,  1.0,  2.0  and  3.0  ksi.  Strains  were 
obtained  only  at  0,  1.0  and  2.0  ksi.  The  compressive  stress  (difference) 
versus  axial  shortening  and  transverse  extension  strains  are  shown  in 
Figure  15,  with  each  curve  the  average  of  4-6  tests  and  scatter  as  shown 
by  bars.  Increased  strength  is  shown  with  confining  pressure,  quantitatively 
agreeing  with  other  data  available  (Ref.  5).  At  2  ksi  confining  pressure 
ductile  behavior  is  shown,  and  relatively  largt  compressive  strains  are 
possible.  The  negative  ratio  of  transverse  to  axial  strains  (the  apparent 
Poisson  ratio)  and  the  slope  of  the  stress  difference-shortening  strain  curve 
(the  apparent  Young's  elastic  modulus)  near  zero  stress  difference  increase 
slightly  with  confining  pressure. 

The  data  from  the  constant  confining  pressure  tests  are  replotted  in 
the  shear  plane  in  Figure  16  and  in  the  dilatation  plane  in  Figure  17. 

In  Figure  16  the  square  root  of  the  second  deviatoric  invariants  of  stress 
and  strain  are  plotted  (equal  to  2//3  times  shear  stress  and  1//3  times 
engineering  shear  strain,  respectively),  with  each  defined  as  shown  below 
in  terms  of  principal  stresses,  o  ,  o  ,  o.  and  principal  strains,  e  , 
c2 ,  c .  (Ref.  33). 

vOj  =  (1/  6)  /(oj*( j2)2  +  (o2-a3)z  +  (a3-a3)2  (1) 

/lj  =  (1//6)  /( G1"e2 )Z  +  (c2~e3^2  +  (Ve3>S  (2) 

The  slope  of  this  curve  is  equal  to  twice  the  apparent  shear  modulus,  G. 

The  values  taken  near  zero  shear  stress  (/J^  <  0.005  ksi)  increase  only 
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Principal  stress-strain  curves  for  compressive 
loading  under  constant  confining  pressure. 
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Figure  16.  Deviatoric  stress-strain  curves  for 
constant  confining  pressure  tests. 
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mean  normal  stress 


Figure  17.  Dilatation  stress-strain  curves  for  constant 
confining  pressure  tests. 
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si ightly  with  pressure. 

In  the  dilatational  plane,  Figure  17,  the  mean  normal  stress,  p  ,  is 
plotted  versus  volume  change,  aV/VqI  defined  as 

pm  =  1/3  ^0l  +  °2  +  c3)  (3) 

aV/Vq  =  Ej  +  e2  +  c5  (4) 

The  slope  of  this  curve  is  equal  to  the  apparent  bulk  modulus,  B.  Pure 
hydrostatic  loading  shows  that  the  hydrostat  tends  to  stiffen  at  very  low 
pressures,  up  to  2-3  ksi,  then  the  concrete  structure  begins  to  collapse 
and  the  hydrostat  softens  until  about  25  ksi  where  it  again  stiffens 
(as  shown  by  the  insert).  Several  pure  hydrostatic  load-unload  tests,  which 
were  conducted  to  determine  at  what  pressure  permanent  compaction  begins, 
indicated  that  measurable  permanent  compaction  begins  at  about  6  ksi  pressure. 

Unconfined  and  constant  confining  pressure  tests  show  dilation  begins 
at  some  load  below  maximum  stress;  however,  before  large  dilation  begins, 
compaction  in  excess  of  the  hydrostatic  compaction  occurs.  This  compaction 
is  believed  due  to  collapsing  of  the  structure  caused  by  the  shear  stress, 
and  has  been  observed  on  other  porous  materials  (Ref.  28). 

In  tension,  bulking  is  observed,  with  the  volume  change  increasing  toward 
greater  volume  increase  with  increased  load.  To  simulate  the  bulking  in 
tension  and  dilation  in  compression  the  apparent  Poisson  ratio  would  have  to 
vary  from  negative  to  greater  than  0.5.  These  curves  show  that  the  assumption 
of  a  linear  elastic  material  with  a  fixed  Poisson  ratio  is  not  reasonable 
even  for  the  low  stresses  shown  here. 

Variable  Load  Path  Tests:  In  addition  to  the  unconfined  and  constant  confining 
pressure  tests,  other  load  path  tests  were  performed,  including  extension, 
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constant  shear  stress  (proportional  to  /d'2) ,  constant  mean  stress  (proportional 
to  Ji),  and  uniaxial  strain. 

The  load  paths  for  these  tests  are  shown  in  Figure  18,  where  stress 
difference,  aa-ot. is  plotted  versus  the  transverse  stress  ot.  Each  curve 
represents  a  single  test.  For  the  constant  Jlt  /J'2  and  uniaxial  strain  tests, 
all  three  stresses  are  compression,  while  for  the  extension  tests  the 
transverse  stresses  are  compression  while  the  axial  stress  may  be  compression 
or  tension. 

Figure  19  shows  the  stress  difference  plotted  versus  axial  shortening 
and  transverse  extension  strains;  that  is,  the  strains  during  the  hydrostatic 
loading  part  of  the  test  have  not  been  included.  The  strains  during  the 
hydrostatic  loading  can  be  obtained  from  Figure  17  by  assuming  axial  and 
transverse  strains  to  be  equal  and  thus  given  by  1/3  of  the  volume  change. 

Figure  20  shows  the  tests  plotted  in  the  mean  normal  stress-volume  change 
plane,  and  Figure  21  shows  the  tests  plotted  in  the  deviacoric  invariant 
str  'ss-strai n  plane. 

Fa i 1 ure  Envelope :  The  magnitude  of  the  deviatoric  stress  at  the  brittle 
maximum  stress  or  at  "he  ductile  yield,  which  is  defined  as  the  stress  at 
]%  axial  strain,  is  plotted  versus  mean  normal  stress  in  Figure  22.  A 
variety  of  load  paths  are  shown,  including  triaxial  compression  and  extension, 
constant  mean  normal  stress,  constant  deviatoric  invariant  stress  and  uniaxial 
strain.  The  triaxial  compression  tests  form  an  upper  bound  while  the  two 
extension  tests  suggest  that  triaxial  extension  forms  a  lower  bound.  The 
uniaxial  strain  load  path  yields  at  a  point,  between  the  triaxial  compression 
and  extension  bounds.  The  'failure  envelope'  suggests  that  for  selected  loads, 
such  as  triaxial  compression,  deviatoric  invariant  yield  stress  (or  yield 
shear  stress)  may  be  related  to  the  first  stress  invariant  (mean  normal  stress). 
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Figure  19.  Principal  stress-strain  curves  for  various  load  paths. 
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Figure  21.  Deviatoric  stress-strain  curves  for  various  load  paths 
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Maximum  stress  locus. 


but  that  for  more  general  loading,  a  unique  relationship  does  not  exist. 
Unloading  Response:  Thus  far,  the  test  data  show  loading  only,  with  tests 
generally  conducted  until  strain  gages  failed  or  (rapid)  fracture  occurred. 

A  number  of  tests  were  specifically  conducted  to  define  the  unloading 
response,  after  some  amount  of  inelastic  deformation  hod  occurred.  Figures 
23-26  show  stress-strain  curves  for  a  number  of  load-unload  paths.  The  data 
are  plotted  showing  mean  normal  stress,  (l/3)(c1+o2+o3) ,  versus  volume 
change,  c]+c2+e 3;  as  square  root  of  the  second  deviatoric  invariant  of  stress 
versus  strain;  and  as  stress  difference,  |oi-oj|,  .ersus  strains,  and  e3. 
Here  oj,  t\  refer  to  axial  direction  and  o2=0 3  and  £3  refer  to  transverse 
direction. 

The  unloading  curves  show  that  unloading  paths  do  not  follow  an  'elastic' 
unload  path,  but  instead  suggest  both  an  anisotropic  hardening  that  causes 
yielding  during  unloading  and  strain  induced  anisotropy  of  the  elastic 
constants.  The  former  is  suggested  by  the  complex  unload  paths  as  the  axial 
compression  is  removed  at  constant  confining  pressure,  while  the  latter  Is 
suggested  by  the  shape  of  the  unload  'hydrostat'  as  the  confining  pressure 
is  removed. 

Reloading  follows  nearly  the  previous  unload  path  If  the  unload-reload 
occurs  under  nearly  the  same  load  ratio.  If  reload  occurs  for  a  grossly 
different  load  ratio,  as  for  example  compression  under  constant  confining 
pressure  followed  by  compression  without  confining  pressure,  then  the  reload 
path  does  not  follow  the  previous  unload  path. 

Fracture:  Fracture  patterns  are  shown  in  the  photographs  of  Figure  27.  The 
photographs  represent  specimens  taken  either  to  brittle  fracture  (well  beyond 
maximum  stress)  or  to  large  ductile  strains,  which  generally  produced  macro 
shear  failures.  Figure  27a  shows  specimens  which  viere  subjected  to  compressive 
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Figure  23.  Triaxial  compression  load-unload  at  constant  confining 
pressure  of  about  2  ksl. 
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loads.  The  unconfined  test  shows  brittle  behavior,  while  1.0  ksi  confining 
pressure  and  above  clearly  show  ductile  response.  The  constant  mean  normal 
stress  (constant  Jj)  and  the  constant  deviatorlc  invariant  stress  (constant 
/j'2)  tests  were  at  sufficiently  high  enough  confining  pressure  so  that  ductile 
behavior  occurs.  Under  uniaxial  strain  loading,  ductile  yield  is  observed. 

Figure  27b  shows  specimens  which  were  subjected  to  tension  or  triaxial 
extension  stress.  A  brittle  behavior  is  shown,  even  for  an  extension  test  at 
10  ksi  confining  pressure;  however,  the  inserts  k  and  1  show  the  different 
failure  mode  between  tension,  where  failure  occurs  almost  totally  due  to 
fracture  in  the  cement  matrix  and  aggregate-cement  separation,  and  extension 
at  high  confining  pressure,  where  failure  occurs  due  to  fractures  of  the 
cement  matrix  and  fracture  of  the  aggregate.  For  both  tension  and  extension 
tests,  fracture  occurred  in  the  specimen  mid-section,  suggesting  that  the 
bonded  end  tabs  were  not  producing  an  unfavorable  stress  concentration  at 
the  ends. 
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SECTION  V 


THEORETICAL  FORMULATION  OF  MODEL 

The  features  of  the  concrete  stress-strain  behavior  observed  in  the 
laboratory  data,  particularly  permanent  set  and  hysteresis  on  unloading, 
suggest  that  a  mathematical  representation  can  be  based  on  a  plasticity 
model.  This  has  been  accomplished  in  that  a  representation  of  the  stress- 
strain  behavior  has  been  developed  that  fits  the  data  reasonably  well. 

The  plasticity  model  employed  uses  a  strain-hardening  yield  surface,  with 
two  segments,  and  an  associated  flow  rule. 

The  general  characteristics  of  the  concrete  mechanical  behavior  are 
a  strong  dependence  of  stress-strain  response  and  maximum  stress  (or  fracture) 
on  mean  normal  stress,  nonlinearity,  hysteresis,  permanent  set  on  unloading, 
and  a  coupling  of  the  shear  and  dilatational  behavior.  Similar  mechanical 
behavior  has  been  observed  in  rocks  and  mathematical  representations  have 
been  studied  previously  that  were  based  on  plasticity  models  (Refs.  23, 

25,  34).  The  volume  expansion  due  to  shearing  seen  in  much  of  the  concrete 
data  is  a  well-known  feature  for  many  materials  a. id  has  been  often  termed 
dilatancy  in  the  soil  mechanics  and  rock  mechanics  literature.  The  behavior 
seen  in  the  mean  stress  versus  volume  strain  curves  for  concrete  is  complex, 
exhibiting  both  compaction  and  bulking  due  to  shearing  in  different  stress 
ranges  in  the  same  test.  This  can  easily  be  seen  in  the  comparison  illustrated 
in  Figure  28.  In  this  figure  the  mean  stress  versus  volume  strain  curve 
of  a  test  without  shearing  (a  hydrostatic  compression  loading)  is  compared 
with  the  same  curve  from  a  triaxial  compression  test.  As  illustrated  in 
the  figure,  a  compaction  first  occurs  due  to  the  shearing  stress  followed 
by  bulking  or  volume  expansion.  This  same  behavior  has  been  observed  recently 
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Figure  28.  Illustration  of  compaction  and  bulking  seen  in 
volume  beha/ior  of  porous  materials. 


in  tests  of  porous  rocks  (Refs.  28,35).  In  addition,  a  relatively  large 
amount  of  compaction  can  be  observed  In  hydrostatic  pressure  tests  at  higher 
mean  ‘.v*«s  levels. 

•  following  discussion  a  plasticity  model  used  to  represent  this 
mechanical  behavior  will  be  described.  The  mooel  uses  an  associated  flow 
law  and  is  based  on  a  strain  hardening  yield  surface  with  two  separate 
segments.  One  of  the  surfaces  is  in  the  form  of  a  "cap,"  as  used  previously 
by  Sandler,  et  al.  (Ref.  24).  The  reasons  for  the  selection  of  the  yield 
surface  in  this  •form  will  be  discussed  below. 

Multiple  Yield-Surface  Model:  In  the  theory  of  plasticity  it  is  assumed 
that  material  behavior  is  elastic  until  the  applied  stresses  reach  a  critical 
value.  Continued  loading  result:  in  both  elastic  and  Irrecoverable  plastic 
straining.  The  critical  value  of  stress  may  be  considered  to  be  a  surface 
in  stress  space  and  thus  is  referred  to  as  a  yield  surface. 

The  yield  surface  to  be  employed  for  concrete  is  illustrated  in  Figure  29. 
Two  independent  segments  are  employed,  termed  Fj  or  __  as  shown.  The  use 
of  two  segments  appears  advantageous  for  modeling  the  volume  plastic  strain. 

As  discussed  previously,  the  concrete  can  show  either  compaction  or  bulking 
in  the  plastic  ( irrecoverable)  volume  strain.  It  can  be  shown  that  a  yield 
surface  with  a  negative  slope  as  the  segment  Fj  (note  that  the  abscissa  of 
Figure  29  is  negative  tu  the  right)  will  snow  plastic  volume  bulking.  Con¬ 
versely  a  .yield  surface  with  the  opposite  slope,  such  as  Fp,  will  shew  plastic 
volume  compaction.  Thus  the  bulking  behavior  is  governed  by  the  segment  Fj 
and  the  compaction  by  F? . 

The  theory  of  multiply  segmented  yield  surfaces  has  been  developed  by 
Koi ter  (Ref.  36)  and  may  be  considered  a  part  of  the  classical  theory  of 
plasticity.  It  should  be  noted  that  the  overall  yield  surface  does  not  have 
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Figure  29.  Illustration  of  two  yield  surfaces  used  for  porous 
material  pi astici ty  model . 


bb 


to  be  smooth.  For  example,  as  Illustrated  in  Figure  29  the  segments  Fi 
and  F2  are  not  tangent  to  each  other  at  their  Intersection. 

A  number  of  assumptions  necessarily  have  to  be  made  to  carry  out  the 
formulation  of  the  model.  The  use  of  the  stress  invariants  /Tj  and  Jj  to 
describe  the  yield  surface  is  a  simplifying  assumption  that  is  inaccurate 
in  fine  detail.  This  assumption  will  be  discussed  further  In  Section  VI. 

It  Is  also  assumed  that  tension  stress  states  can  be  handled  similarly  to 
compression,  with  the  decrease  in  /Jj  as  illustrated  in  Figure  29 
accounting  for  the  weakness  in  tension.  Further  assumptions  are  that  an 
associated  flow  rule  and  isotropic  hardening  can  be  used  in  the  plasticity 
model.  The  use  of  an  associated  flow  rule  relates  the  plastic  strain 
Increments  to  the  shape  of  the  yield  surface  in  a  specific  manner  which 
will  be  made  explicit  below.  It  is  advantageous  in  the  sense  that  it 
greatly  simplifies  certain  theoretical  problems  such  as  whether  or  not 
the  stress-strain  law  is  thermodynamically  admissable.  Isotropic  hardening 
means  that  the  yield  surface  expands  uniformly  as  the  material  strain 
hardens.  This  assumption  in  effect  Insures  that  the  unloading  behavior 
will  always  be  elastic,  which  apparently  is  not  appropriate  In  detail  for 
concrete.  However  isotropic  hardening  is  at  least  a  first  step  and  can  be 
considered  an  approximation  to  the  real  unloading  behavior.  In  general 
these  assumptions  can  be  changed  as  detailed  investigations  of  the  various 
aspects  of  mechanical  behavior  are  carried  out. 

In  the  theory  of  plasticity  for  small  strains  the  strain  increment  can 
be  separated  into  elastic  and  plastic  components  as 


The  elastic  strain  increment  component  is  related  to  the  stress  increment 
through  the  elastic  stress-strain  law.  The  plastic  strain  increment  can 
be  expressed  by  the  associated  flow  law  as  (Ref.  37) 

de.-.P  =  sF  a  (6) 

*  J 

3oij 

where  f  is  the  yield  function  and  a  is  a  constant  that  for  strain  hardening 
materials  will  contain  the  stress  increments,  thus  providing  a  relationship 
between  the  plastic  strain  increments  and  the  stress  increments. 

Koiter  (Ref.  36)  has  shown  that  the  associated  flow  rule  for  a  yield 
surface  with  two  segments  can  be  written  as 


de^jP  a  3 F j  Aj  +  3F?  A?. 


(7) 


where  Fj  and  F2  are  the  two  yield  surface  segments.  It  can  be  seen  that 
equation  (7)  is  similar  to  (6)  except  that  plastic  strain  increments  can 
be  obtained  from  both  yield  surface  segments.  It  will  be  shown  later  that 
Aj  and  A 2  will  both  be  nonzero  only  if  the  state  of  stress  is  that  defining 
the  intersection  of  the  two  yield  segments.  At  this  point  the  plastic  strain 
increment  is  the  sum  of  contributions  from  both  yield  surface  segments.  This 
explains  why  the  intersection  does  not  have  to  be  smooth.  It  is  only 
necessary  that  each  segment  be  individually  smooth  so  that  the  terms 


dFj  and  aF? 

3oi  j  Sc,1j 

each  exist. 

The  constants  Aj  and  can  be  defined  by  specifying  appropriate  strain 
hardening  rules.  This  will  be  carried  out  in  general  form  here  and  the 
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specific  forms  necessary  to  fit  the  experimental  data  will  be  examined  in 
Section  VI. 

The  yield  surface  segments  are  assumed  to  have  the  following  forms 


=  ^  +  9i  (Jj)  “  Kj 

(8) 

(^•)2  +  r2  (J,)*  -  K2 
~2~  1 

(9) 

where  gt  (JJ  is  a  function  to  be  determined,  and  r  is  a  constant.  As  usual, 
ZZ'  is  the  second  stress  deviator  invariant  and  is  the  first  stress 
invariant  given  by 


Jj  a  Jj  S 


ij  1J 


33)2+(a33~all>2  +  a122  +  °2  32  +  °312 

6 


and 


J,  =  o  ,  .  +  n  +  n 

A  1  i  ZZ  J  J 


The  initial  values  of  the  constants  «i  and  K2  define  the  locations  of 
the  yield  surface  segments.  When  the  stress  state  is  located  inside  the 
yield  surface  the  stress-strain  response  is  elastic,  thus  Fj  =  0  or  F2  =  0 
is  a  necessary  condition  for  yielding. 

In  a  strain  hardening  material  the  yield  surface  can  move  as  plastic 
straining  occurs.  This  can  be  accounted  for  by  changing  the  values  of  the 
constants  Kj  and  «2  as  plastic  straining  takes  place.  The  relationship  of 
the  change  in  the  constants  to  the  straining  is  termed  the  strain  hardening 
rule.  Since  in  the  present  case  it  is  desired  to  have  the  two  yield  surface 
segments  move  independently  it  is  necessary  to  have  a  hardening  rule  for 
each  segment.  Once  the  form  of  the  yield  function  is  established  the  plastic 
stress  strain  response  is  fully  determined  by  the  hardening  rule  employed. 
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It  was  assumed  that  the  movement  of  the  segment  Fj  depended  on  the 
plastic  shearing  strain  and  that  the  movement  of  F2  depended  on  the  plastic 
compaction  strain.  Specifically  the  hardening  rules  were  taken  as 


d/T'P  =  hj  dKj 

(10) 

devc^  =  ^2  ^2 

(ID 

where  d /T'2^  is  the  second  invariant  of  the  plastic  strain  increment  tensor, 
drvcp  is  the  volume  strain  increment  associated  with  compaction  only,  and  hj 
and  h2  are  functions  to  be  determined  by  the  experimental  data.  The  total 
plastic  volume  strain  has  thus  been  divided  into  two  parts  as 

dcvp  -  (dtV{j  ) b ul k i ng  +  ^cvc^compaction 

This  distinction  in  the  plastic  volume  strain  is  made  for  purposes  of  the 
derivation  only  and  will  not  appear  in  the  final  stress  strain  law. 

The  assumption  of  the  form  of  the  hardening  rules  (equations  10  and  11) 
is  difficult  to  justify  until  more  is  known  about  microstructural  response 
mechanisms  of  concrete.  The  use  of  the  deviatoric  invariant  in  equation  10  is 
standard  in  the  plasticity  literature  for  metals  (Refs.  33,37)  and  has  been 
used  for  geologic  materials  (Refs.  4,23).  The  use  of  the  compaction  strain  in 
equation  11  appears  desirable  since  the  yield  segment  F2  should  harden  (expand) 
as  compaction  occurs  In  both  hydrostatic  and  nonhydrostatic  loadings.  However, 
it  should  be  recognized  that  these  are  plausibility  arguments.  While  the 
appropriateness  of  the  assumptions  of  equations  10  and  11  can  in  part  be 
examined  by  the  fit  between  the  resulting  model  and  experimental  data,  the 
full  validity  oi  extrapolation  to  new  stress  conditions  is  not  establisned. 
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The  constants  Aj  and  A?  are  determined  by  substituting  the  hardering 
rules  into  the  appropriate  equations  as  follows.  Since  plastic  shearing 
strains  are  obtained  from  loading  on  both  yield  segments,  the  plastic  strain 
obtained  from  both  segments  (equation  7)  should  be  substituted  in  equation  in. 
This  substitution  is  straightforward  but  lengthy  and  is  given  in  Appendix  C. 
The  result  is 

3 F j  Aj  +  3F2  A 2  =  2h]  d K i  (13) 

d  </j 2  3  SO*, 

The  compaction  plastic  volume  strain  is  obtained  only  from  the  F2 
yield  surface  segment,  provided  that  the  slopes  of  the  segments  Fj  and  F2 
ar>_  restricted  appropriately.  With  this  implicit  restriction,  the  strain 
from  the  F2  segment  which  is  given  as 

(d£i,p)  -  aF2  A2  (14) 

30ij 

should  be  substituted  into  the  hardening  rule  for  compaction,  equation  11. 

The  result  of  this  substitution  gives 

3F2  a2  =  1/3  h2  dK2  (15) 

fjJ  j 

To  ensure  that  plastic  straining  occurs  only  when  the  material  is 
yielding,  it  is  necessary  to  define  the  constants  A,  and  a2  as  follows: 
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■  ,■  at :I  n::n"|;  | ,  ''  . '  i '  T:p"l r  'hi  \\>W  i|7'T",i:n,  jifia?  ft  HT  ^ 1 


Ai  if  Fx  =  0 
dF)>0 

0  if  Fx<0 

or  dF^O 


' 

A 2  if  ?2  ~  0 
\d?2>0 

x2  = 

0  if  ( F?<0  ,  % 

(or  dF2$0  (17) 

Equation  15  can  then  be  solved  for  a2,  subject  to  equation  17  above,  as 

1/3  h2  dK2 

A  2  =  ^2  (18) 

3  J ) 

and  equation  13  can  be  solved  as 


_J__  (2hl  dKj 
3F~i 

3/T2 


3F2  A2) 

d75'z 


It  should  be  noted  that  the  yield  surface  segments  do  not  move  inward 
during  any  loading.  This  observation  follows  directly  from  the  hardening 
rules  assumed  in  equations  10  and  11,  and  the  assumption  of  isotropic 
hardening.  It  is  assumed  that  the  plastic  deviatoric  strain  invariant  and 
the  plastic  compaction  strain  are  non-decreasing  in  absolute  magnitude. 

Thus  the  functions  hi  and  h2  are  restricted  to  not  changing  sign  during  yielding. 

The  above  equations  comolete  the  basic  derivation  of  the  two  yield 
surface  segment  plasticity  model.  In  Section  VI  the  necessary  functions  and 
constants  will  be  fitted  to  experimental  data  on  concrete  and  the  basic  utility 
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of  the  model  will  be  demonstrated. 


Derivation  of  Strain  Increment  as  Independent  Variable:  The  usual  plasticity 
model  formulation  as  shown  above  considers  the  strain  increments  resulting 
from  a  stress  increment  loading,  so  that  the  stress  increments  can  be  considered 
to  be  the  independent  variables.  In  continuum  analysis  codes,  particularly 
those  using  the  finite-element  formulation,  it  is  often  necessary  to  have 
the  incremental  stress-strain  relationship  in  the  inverse  form  so  that  the 
strain  increments  are  the  independent  variables.  This  reformulation  can  be 
readily  carried  out  without  introducing  any  changes  in  the  assumptions  of 
the  material  behavior.  This  reformulation  is  carried  out  as  follows: 

Consider  the  usual  separation  of  strain  increments  into  elastic  and 
plastic  components  stated  previously  as 


+  dcP . 
ij 


The  elastic  stress-strain  law  can  be  written  as 


doij =  Ld?kk  6ij +  2Gdc'j 

where  the  symbol  L  is  used  here  for  the  Lame  constant 


(5) 


(20) 


L  s  v  E 

(l+v)(l-2v) 


In  inCex  notation  this  can  be  written  as 

doi j  =  Ci jkl  dEkl  (21* 

where  C  is  the  elastic  coefficient  matrix.  Substituting  equation  5  into 
21  gives 


do.  .  =  C. .. ,  dc.  .  -  C. 

ij  ijkl  kl  i jkl  kl 


(22) 
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Using  the  flow  rule  developed  for  the  two  segment  yield  surface 


de[!]  =  3  F.  Ai  +  3  F2  A2 


then  equation  22  can  be  written  as 


doi j  =  Ci jkl  dckl  "  Ci jkl  '  3Fl  Xl  +  3F?  *2 


The  hardening  rules  were  expressed  previously  as 


d/l'2P  =  h,  dFj  =  hi  3 F x  do  — 


d £yC  =  h2  dF2  =  h2  dF 2  do.j 


where  dePc  is  the  volumetric  compaction  from  the  F2  yield  surface.  The 
equations  previously  developed  for  Aj  and  a2  were 


DF,  a,  dF2  A2  =  2  h,  3 F j  da-j j 


3F2  a 2  =  1/3  h2  dF2  do^j 


Substituting  equation  23  for  da.j  into  the  right  hand  sides  of  equations 
13  and  15  above  gives 


dF,  A,  +  jP.  > =  2  h .  a  F ,  iC.,,dr,,*C 


..  ,  /  aF  ,  A,  +  D  F  A  -A 


’||  ||-1|-^  -,*knil  mi M  ■-IL.Ii.I.lJ - iriit. jJt lW  .llilkl'ull,  1* 
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If  either  Ax  =  0  or  A2  ~  0  because  loading  does  not  occur  on  both  yield 
surfaces,  equations  28  and  29  should  be  used.  Thus,  if  either  <  0 

or  dF2  $  0,  then 


x2  =  0 

xl  “  A3  =  a6 

Aj  A4 

Likewise,  if  either  Fj  <  0  or  dFj  <  0,  then 

Aj  *  0 

A2  =  A3  =  A6 
A?  As 


The  final  stress-strain  law  is  given  in  equation  23  as 


do1  j  =  Cijkl  d£kl  '  Cijkl  / £i_  +  ^ 

\3ak1  dakl  / 


where  x1  and  a2  are  given  as  equations  30,  31  or  32,  33  above. 


SECTION  VI 

COMPARISON  OF  MODEL  AND  EXPERIMENT 

In  this  section  the  functions  in  the  previous  equations  will  be  specified 
to  fit  the  experimental  data  on  concrete.  Stress  is  expressed  in  units  of 
psi  with  the  convention  of  tensile  stresses  taken  as  positive. 

The  elastic  behavior  of  the  concrete  was  taken  as  linear  in  both  shear 
and  bulk  with  constants 

G  =  2.0  x  106  psi  shear  modulus 

B  =  3.0  x  106  psi  bulk  modulus  (34) 

The  yield  functions  were  taken  as 

Fj  =  /J'2  -  1000.  [l2. 2  -  11  exp  (Jj/40000.)]  [l.  -  exp((Jr800.  )/700.  S]-^ 

(35) 

F 2  -  (/J. ‘ V  +  *  J,2  -  K2  (36) 

2 

The  yield  surface  segment  Fj  is  taken  as  the  locus  of  points  where 
deviation  from  the  initial  elastic  and  approximately  linear  behavior  in  shear 
occurs.  This  definition  is  somewhat  arbitrary  as  the  point  of  deviation 
from  linearity  depends  on  subjective  judgment  and  experimental  accuracy. 

However,  the  problem  is  not  serious  if  the  hardening  rule  parameters  are 
specified  accordingly. 

The  segment  Fj  is  based  on  the  tension,  unconfined  compression,  and 
triaxial  compression  test  results.  As  a  smooth  transition  from  tension  to 
compression  is  assumed,  no  specific  provision  for  a  tension  cutoff  is  necessary. 

The  yield  surface  segment  F  is  defined  in  at;  arbitrary  manner.  That  is, 
it  was  taken  as  being  small  with  respect  to  tension  and  compression  failure 
stresses  and  thus  not  well  defined  by  the  experimental  data.  The  hardening 
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rule  was  adjusted,  however,  so  that  the  combined  effect  of  yield  surface 
and  hardening  rule  did  match  experimental  data.  A  comparison  of  equation  36 
with  equation  9  shows  that  the  ratio  of  the  radii  of  the  F2  ellipse  has 
been  taken  as  2.  Trial  runs  showed  that  the  resulting  fit  was  not  sensi¬ 
tive  to  this  parameter,  thus  the  value  of  2  is  somewhat  arbitrary.  The 
constants  Ki  and  K2  change  as  strain  hardening  occurs  according  to  strain 
hardening  rules  given  below  and  have  initial  values 

Kj  =  0.0  psi 

K2  =  5  x  107  (psi)2  (37) 

The  hardening  rules  are  given  as 

d/I]P  =  h1  dFj  =  hj  dK3  if  dFj  >  0  (38) 

with 

hj  =  10~6  J  1  -  100  [1  -  exp(120  /r^p/(l  -10'4  Jj))]  j  (39) 

2  -  7  x  10'5  J, 

and 

d%cP  “  h2  &?.  =  h2  dK2  if  df2  >  0  (40)  . 

wi  th 

h2  =  |  -2.4  HP4  /u2  [  i  -  exp  (-vT2/4000  )]  [exp(-^4/(8  x  1035))J 

-  0.702  x  10" 1 1  [  1  -exp  ( J  j 3/  ( 1 . 35x10 1 3 ) )]  [exp(Jj/(6xl04))]}  (41) 

As  discussed  in  the  previous  section  the  constants  Kj  and  K2  defining  the 
yield  surface  are  changed  as  strain  hardening  occurs  and  are  defined  so  as  to 
be  nondecreasing.  Thus,  each  yield  surface  segment  moves  away  from  the 
origin  as  it  is  loaded.  During  loading  on  F,  (i.e.,  Fj  =  0  and  dF:  0)  dF}  = 
dKj  and  during  loading  on  F2  (i.e.,F2  =  0  and  dF~  >  0)  dF2  =  dK2.  Otherwise, 


dK]  -  dK2  =  0.  Unloading  from  and  reloading  back  to  the  yield  surface  is 
thus  assumed  to  be  elastic.  Changes  in  stress  from  tension  to  compression 
are  defined  by  the  above  procedure,  i.e.,  no  special  or  additional  assumptions 
need  be  made. 

The  stress-strain  behavior  based  on  these  equations  has  been  computed 
and  compared  with  the  laboratory  data  in  Figures  30  through  41.  Figure  30 
shows  a  hydrostatic  loading  comparison.  In  this  figure  the  nonlinearity  of 
the  curve  above  a  pressure  of  about  6,000  psi  is  caused  by  compaction. 

A  comparison  with  unconfined  tension,  unconfined  compression,  and  triaxial 
compression  is  shown  in  Figures  31  and  32.  The  shearing  curves  in  Figure  31 
demonstrate  the  nonlinearity  due  to  yielding.  The  dilatational  response 
clearly  shows  compaction  and  bulking  in  both  the  model  and  experimental  data. 

Other  comparisons  are  shown  in  Figures  33  through  41.  Figures  33  and  34 
show  the  comparisons  for  the  tests  in  which  /I*  and  are  held  constant  over 
part  of  the  loading.  Figures  35  and  36  show  comparisons  for  an  extension 
test  at  a  confining  pressure  of  approximately  5  ksi.  Figures  37,  38,  and  39 
give  the  one-dimensional  strain  test  comparison.  Figures  40  and  41  show 
load  and  unload  triaxial  compression  test  results. 

In  general,  the  comparison  between  model  and  experiment  is  reasonably 
good  and  the  qualitative  features  are  well  represented  in  the  model.  The 
one-dimensional  strain  test  comparisons  are  somewhat  poorer  than  the  other 
tests.  This  can  largely  be  ascribed  to  the  high  mean  stress  levels  in  this 
test.  The  mean  stress  levels  in  the  triaxial  compression  tests  were  lower 
and  the  model  parameters  were  based  primarily  on  these  tests.  Thus,  the 
model  had  to  be  extrapolated  significantly  for  the  one-dimensional  strain 
test.  It  should  also  be  noted  that  the  variations  in  the  experimental  data 
in  the  one-dimensional  strain  test  were  significant.  For  this  reason  the 
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concrete 


MEAN  NORMAL  STRESS  KSf 


Figure 


32.  Comparison  of  model  and  experimental  dllatational  stress-strain 
response  in  constant  confining  pressure  tests. 
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DEV'A^ORiC  INVARIANT  STRESS  KS* 


Figure  3j.  Coinpa ri son  of  model  and  experimental  shearing  stress-strain 
response  in  constant  >rCy  and  constant  J,  tests. 
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Figure  34.  Comparison  of  model  and  experimental  dilatational  stress-strain 
response  in  constant  /Ji,  and  constant  J ]  tests. 


DEV  ATOR'C  INVARIANT  STRESS  -43|  KSl 


'jure  35. 


Comparison  of  model  and  experimental  shearing  stress 
strain  response  in  extension  test. 
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Figure  36.  Comparison  of  model  and  experimental  dilatations!  stress-str 
response  in  extension  test. 


ain 
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Figure  38.  Comparison  of  model  and  experimental  dilatational  stress-strain  response  in 
one-dimensional  strain  tests. 


DEViATORIC  INVARIANT  STRESS  -4j|  KSI 


Figure  39.  Comparison  of  node!  and  experimental  stress  path 
response  in  one-dimensional  strain  tests. 
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TRIAXIAL  STRESS  LOAD/UNLOAD 
TEST  AT  2  KSI  PRESSURE 


Figure  40.  Comparison  of  model  and  experimental  shearing  stress 
strain  response  in  constant  confining  pressure  load/ 
unload  test. 


TRIAXIAL  STRESS  LOACVUNLOAD 
TEST  AT  2  KSI  PRESSURE 


Figure  41,  Comparison  of  model  and  experimental  dilatational  stre 
strain  response  in  constant  confining  pressure  load/ 
unload  test. 


experimental  data  from  both  tests  performed  are  shown  individually  in 
Figures  37  through  39. 

Computer  Subroutine:  The  equations  of  thp  previous  section  have  been  programmed 

as  a  subroutine  that  calculates  a  stress  increment  for  a  given  strain  increment, 

state  of  stress,  and  plastic  strain.  A  list  of  definitions  of  the  subroutine 

variables  and  a  listing  of  the  subroutine  are  given  in  Appendix  A.  In  general 

the  variable  names  have  been  chosen  tn  conform  to  the  usage  of  the  SLAM 

computer  code.  The  subroutine  utilizes  an  index  notation  for  stress  and  strain 

according  to  the  following  notation,  using  stress  as  an  example: 

SIGI  (1)  5  radial  stress,  orr 

SIGI  (2)  s  hoop  stress,  j 

SIGI  (3)  c  axial  stress,  azz 

SIGI  (4)  c  shear  stress,  a 

A  description  of  the  subroutine  as  presently  written  can  best  be 
sunmarized  by  the  subroutine  statement 

SUBROUTINE  PLASTK  (DEPST,  SIGI,  SI2P,  EPSPI,  ZK1,  ZK2) 

The  variables  are  defined  as 

DEPST  (I)  i  Strain  increments 

SI2P  s  Plastic  strain  deviatoric  invariant,  /T'P 

EFSPI  (I)  P’.a Stic  strains 

ZK1  a  Constant  in  first  yield  surface 

ZK2  =  Constant  in  second  yield  surface 

The  strain  increments  are  calculated  in  the  main  program.  The  remaining 
enables  are  appropriately  incremented  in  the  subroutine  and  mnst  be  saved 
in  storage  in  the  main  program.  The  storage  of  SI2P,  the  plastic  strain 
deviutoric  invariant  could  be  eliminated  as  this  parameter  could  be  calculated 
from  the  plastic  strain. 

Examination  of  Model  Parameters:  Although  a  iarg*j  number  of  constants  are 
used  in  the  curve  fitting  expressions,  only  six  material  functions  are 
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actually  involved.  These  are  the  two  linear  elastic  constants,  the  two 
yield  functions,  and  the  two  hardening  rules.  Although  a  systematic  study 
of  how  these  functions  would  vary  with  concrete  manufacturing  variables  has 
not  been  undertaken,  it  appears  reasonable  to  assume  that  the  functions  would 
vary  in  a  systematic  manner  with  changes  in  the  concrete  material  properties. 
In  the  following  an  estimate  of  how  the  functions  could  be  changed  to 
accommodate  different  concrete  properties  is  given. 

The  minimum  test  data  necessary  to  define  the  six  functions  could  consist 
of  stress-strain  data  from  an  unconfined  compression  test,  a  triaxial 
compression  test  at  2  ksi  confining  pressure,  and  a  hydrostatic  loading  test. 
Further  tests  would,  of  course,  add  confirmation  to  the  results.  The  six 
material  functions  would  then  be  obtained  as  follows: 

Bulk  Modulus  B:  The  slope  of  the  approximately  linear  portion, 
say  up  to  a  pressure  of  4  ksi,  of  the  hydrostatic  test  plotted  as 
pressure  versus  volumetric  strain. 

Shear  Modulus  G:  One  half  times  the  slope  of  the  approximately 
linear  initial  portion  of  the  triaxial  test,  plotted  as  versus 

'''I' 

Yield  Function  Fx:  The  function  Fj  is  given  in  the  form  of 

F,  =  Q  +  {g1(J1)  -  K,} 

The  term  in  brackets  {  }  should  be  multiplied  by  the  ratio  of  the  unconfined 
compressive  strength  to  -7000  psi ,  which  is  the  value  of  the  unconfined 
compressive  strength  of  the  concrete  used  in  the  present  studies,  that  is 
multiplied  by  uc/-7000. 
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Yield  Function  F?:  No  change 

Hardening  Rule  hi:  No  change 

Hardening  Rule  h2:  The  function  h2  is  written  In  two  terms  as 

h2  =  -f(/T£)  -  gtJj) 

The  function  f  should  be  multiplied  by  the  ratio  of  the  strain  termed  "S" 
in  the  P  s  2  ksi  triaxial  test  illustrated  in  Figure  42  to  the  same  strain 
in  the  present  concrete,  taken  as  0.0018.  That  is,  multiply  "f"  by  S/0.0018. 

The  function  g  of  h2  should  be  multiplied  by  the  ratio  of  the  strain 
termed  "T"  in  the  hydrostatic  test  illustr.ced  in  Figure  43  to  the  same 
strain  in  the  present  concrete,  taken  to  be  0.029.  That  is,  multiply  "g" 
by  T/0.G29. 

It  should  be  emphasized  that  the  above  procedure  has  not  been  verified 
experimentally.  Rather  it  represents  an  estimate  of  how  the  model  of  the 
present  investigation  might  be  modified  for  a  somewhat  different  concrete. 
Discrepancies  Between  Model  and  Data:  Certain  features  of  concrete  behavior 
exhibited  in  the  laboratory  have  not  been  expressed  explicitly  in  the 
mathematical  representation  at  this  time.  One  of  these  features  is  a  path 
dependent  unloading  seen  in  the  dilatational  stress-strain  curves  of  tests 
unloaded  from  near  maximum  load.  In  the  model,  elastic  unloading  is  assumed 
in  both  shear  and  bulk.  The  real  material  behavior  Indicates  that  part  of 
the  volumetric  bulking  may  be  recoverable,  and  in  a  way  that  depends  on  the 
unloading  path.  An  example  of  unloading  in  a  triaxial  (constant  confining 
pressure)  test  is  shown  in  Figures  40  and  41.  The  model  assumes  isotropic 
behavior'  and  does  not  include  the  anisotropy  produced  by  directional 
microcracking.  This  behavior  would  show  up  In  reloading  tests  of  concrete 
loaded  near  maximum  stress,  and  particularly  if  the  directions  of  the 
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Figure  42.  Illustration  of  volumetric  strain  "S"  used  for  correction 
of  model  parameter  H2  for  different  concrete. 
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Figure  43.  Illustration  of  volumetric  strain  "T"  used  for  correction 
of  model  parameter  H.,  for  different  concrete. 


principal  stress  axes  were  rotated. 

Another  area  of  uncertainty  1 n  the  present  model  is  in  the  detailed 
description  of  the  yield  surface  and  maximum  stress  surface.  The  maximum 
stress  locus  has  been  represented  in  versus  J1  coordinates  in  the  model. 

It  has  been  well  established  in  the  literature  and  in  the  present  program 
that  the  maximum  stress  locus  is  not  unique  in  these  coordinates,  but  can 
vary  with  the  type  of  test.  In  general  as  discussed  earlier,  the  triaxial 
extension  arid  triaxial  compression  ten  c  results  appear  to  form  separate 
curves,  with  biaxial  compression  results  falling  in  between.  Various 
techniques  have  been  used  to  rectify  this.  Mills  and  Zimmerman  (Ref.  20) 
have  modified  the  shearing  invariant  with  the  addition  of  a  "rotational 
term."  In  a  similar  situation  for  rocks,  Mogi  (Refs.  38,39')  have  modified 
the  mean  normal  stress  invariant  by  adding  a  weighting  factor  to  the  inter¬ 
mediate  principal  stiess.  An  example  of  this  for  granite  is  shown  in  Figure  44. 
It  should  be  noted  that  although  the  fit  is  apparently  quite  consistent 
in  the  coordinates  of  Figure  44,  a  replotting  of  the  data  in  terms  of  a 
biaxial  stress  envelope  can  lead  to  anomalies.  Further  complicating  the 
situation  is  the  normality  rule  (associated  flow  rule)  used  in  the  plasticity 
model  that  serves  to  link  the  strength  and  stress-strain  behavior.  Although 
some  biaxial  stress-strain  data  are  available,  as  by  Kupfer,  Hilsdorf,  and 
Rusch  (Ref.  15)  for  example,  the  situation  has  not  been  adequately  studied. 
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Figure  44.  Failure  locus  for  Westerly  granite  (uniform  stress). 


SECTION  VII 
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COMPARISON  TO  THICK  WALLED  CYLINDER  TESTS 

A  series  of  tests  was  performed  on  thick-walled  concrete  cylinders. 
These  tests  consisted  of  combinations  of  axial  compression  and  internal 
pressure,  producing  combined  compressive  and  tensile  stresses  in  the 
concrete.  To  ensure  an  adequate  wall  thickness  to  aggregate  ratio, 

3/8-inch  nominal  aggregate  size  was  used  with  a  2-inch-thick  cylinder  wall. 
The  overall  cylinder  configuration  was  7  inches  I.D.,  11  inches  O.D. , 
and  20  inches  long. 

Since  the  concrete  constitutive  model  was  formulated  before  the 
cylinder  tests  v/ere  performed,  these  tests  served  as  a  verification  check 
of  model  parameters.  A  description  of  the  analysis  and  comparison  of  theory 
and  experimental  results  follows. 

Analysis ;  If  the  concrete  cylinders  are  treated  as  "thin  walled",  the  hoop 
stresses  are  given  simply  as 


°ea  =  £!1 

where  p  is  the  internal  pressure,  r  is  the  radius  (presumably  an  average 
of  the  inside  and  outside  radius),  and  t  is  the  wall  thickness.  However, 
because  the  radius  varies  from  3.5  inches  to  5.5  inches  it  was  considered 
necessary  to  use  thick-walled  cylinder  theory  to  determine  the  distribution 
of  radial  and  circumferential  stresses  and  strains.  Th i ;  was  accomplished 
by  developing  a  finite  difference  computer  code  that  used  the  concrete 
elastic-plastic  constitutive  model.  The  equations  were  derived  by  considering 
the  usual  equations  (Ref.  AO). 
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equilibrium:  do„  +  a  -a  -  0  (42) 

_!ZL  rr  ee 

dr  r 

strain-displacement:  e  =  dU  (43) 

dr 

£ee  "  U/r 

stress-strain  law  (in  incremental  form); 

do  =  C  de  (44) 

where  C  is  the  concrete  elastic-plastic  stress-strain  matrix.  These 
equations  were  written  in  finite  difference  form  after  eliminating  the 
strains  in  terms  of  the  radial  displacement  U.  The  finite  difference 
equations  were  readily  programmed;  a  listing  of  the  program  is  given  in 
Appendix  R.  The  axial  stress  and  internal  pressure  were  input  in  incremental 
form.  The  axial  stress  was  assumed  uniformly  distributed  throughout  the 
cylinder  and  related  to  the  measured  axial  load  by  the  relationship 

ozz  -  F/A  (45) 

A  check  solution  is  shown  in  Figure  45  for  the  case  of  internal  pressure 
only  using  elastic  properties.  It  can  be  seen  that  good  agreement  was 
obtained  using  11  mesh  points.  This  number  was  used  for  all  of  the  concrete 
cylinder  solutions. 

Discussion  of  Results:  The  results  of  the  cylinder  tests  are  shown  in  Figures 
46  through  52.  In  these  figures  the  axial  stress,  internal  pressure,  axial  and 
circumferential  strain  both  at  the  inner  and  outer  surface  of  the  cylinder  are 
plotted  versus  the  elapsed  time  of  the  test.  Also  shown  in  these  figures  are  the 
strains  in  the  axial  and  circumferential  directions  at  the  inner  and  outer  surface 
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predicted  by  the  numerical  solution. 

In  general  the  comparison  between  the  numerical  solution  and  the 
experimental  results  is  within  :he  consistency  of  the  test  results.  The 
comparisons  for  tests  8*3  and  3-4  shown  in  Figures  51  and  52,  respectively/ 
were  relatively  much  poorer  than  the  others.  It  is  not  apparent  why  this 
should  be;  however,  both  of  these  tests  had  a  very  high  initial  axial  pre¬ 
load  before  the  internal  pressure  was  applied.  Apparently  the  model 
exaggerates  the  circumferential  strain  at  these  stress  levels.  This  could 
indicate  that  the  yield  function  Fj  should  be  increased  somewhat  more 
for  the  cyli'Mer  concrete,  or  else  the  dilatancy  effect  may  be  too  pro¬ 
nounced  for  this  stress  condition.  However,  the  relatively  much  better 
comparison  for  the  other  tests  suggest  that  the  experimental  results  for 
tests  8-3  and  8-4  may  not  be  entirely  consistent  with  the  others. 

An  interesting  phenomenon  was  noted  in  the  analysis  of  the  cylinders. 

In  the  usual  elastic-plastic  cylinder  under  internal  pressure,  using  say  a 
Tresca  or  Von  Mises  yield  criterion,  yielding  starts  at  the  inner  core  and 
progresses  outward  as  the  pressure  is  increased-.  In  the  present  analysis 
yielding  was  observed  to  start  at  or  near  the  outer  surface.  This  can  be 
(...plained  by  considering  the  pressure  sensitivity  of  yielding  in  concrete. 

The  pressure  at  the  inner  bore  was  evidently  sufficient  to  shift  the  location 
of  initial  yielding.  As  a  consequence  of  this,  the  "thick-wall"  effect  is 
not  as  pronounced  for  the  concrete  analysis  as  it  would  be  for  a  pressure- 
independent  elastic-plastic  material. 

Summary:  Axial  force  and  internal  pressure  tests  were  run  on  concrete 
cylinders.  These  cylinders  were  analyzed  using  thick-walled  cylinder  theory 
in  conjunction  with  the  material  property  model  developed  for  concrete. 

Comparison  of  theory  and  experiment  was  relatively  good  for  all  tests, 
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Figure  46.  Comparison  of  model  and 


experiment  for  cylinder  7-1. 
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Figure  52.  Comparison  of  model  and  experiment  for  cylinder  No.  8-4. 


although  cylinders  8-3  and  8-4  were  much  poorer  than  the  others.  The 

evidence  for  explaining  the  poorer  comparison  of  these  two  tests  Is  not 
conclusive. 
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SECTION  VIII 


SUMMARY  AND  CONCLUSIONS 

1.  Although  comparison  of  test  data  from  different  laboratories  is 
very  difficult  because  of  the  large  variability  and  the  complex  response 

of  concrete,  data  available  in  the  literature  do  show  some  of  the  deformation 
features  of  concrete,  including  the  increase  in  strength  with  pressure, 
the  dependence  of  the  yield  on  stress  state  and  dilation  which  occurs  near 
maximum  stress.  Tests  here  provided  added  detail  to  the  stress-strain  response 
and  showed  added  features  of  deformation  including  a)  compaction  associated 
with  shear  deformation  occurs  before  the  onset  of  large  scale  dilation  and 
b)  the  unloading  is  not  elastic  but  shows  complex  unload  paths  which  suggest 
anisotropic  hardening  and  strain  induced  anisotropy.  To  a  lesser  extent 
the  tests  here  showed  the  extreme  sensitivity  of  the  mechanical  properties 
to  casting  procedures,  and  he  importance  of  the  specific  aggregate  surface 
area  in  determining  concrete  strength.  The  strength  appeared  to  be  far 
more  sensitive  to  surface  area  than  to  aggregate  size  variation  from  3/16 
to  3/4 -inch  diameter  aggregate. 

2.  The  r  ndorn  and  relatively  homogeneous  distribution  of  aggregate, 
cement  matrix,  water  and  voids  and  the  general  response  of  the  concrete, 
including  dilation  and  permanent  set,  are  suggestive  of  a  continuum,  elastic- 
plastic  model.  It  is  not  implied  that  concrete  exhibits  'metal  plasticity', 
but  that  plasticity  models  might  adequately  be  used  to  phenomenologically 
model  concrete  as  has  previously  been  done  very  successfully  for  geologic 
materials.  In  a  sense  the  model  used  here  is  an  extension  of  the  initial 
plasticity  model  used  by  Swanson  to  fit  rock  data,  with  the  yield  cap  as 
suggested  by  Baron,  et  al.  An  associated  flow  law  and  a  segmented  yield 
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function  have  been  used  with  a  hardening  law  which  includes  the  shear  and 
dilatation  terms.  The  segmented  yield  and  selected  hardening  laws  allow 
for  compaction  followed  by  dilation,  as  well  as  the  more  conventional  pressure 
sensitivity  of  yield. 

3.  An  Independent  check  of  the  model  has  been  made  by  calculating  surface 
strains  on  a  thick-walled  cylinder  subjected  to  a  programmed  internal  pressure 
and  axial  load,  and  comparing  the  calculated  strains  to  those  measured 

during  the  loading.  The  thick-walled  cylinders  represented  a  concrete  with 
a  somewhat  different  strength,  and  hence,  in  effect  represented  a  different 
batch  of  concrete.  Comparison  between  the  finite  difference  calculation  and 
measured  strain  was  relatively  good  and  provided  some  confidence  r.iiat  the 
model  developed  to  generally  fit  the  features  of  deformation  of  concrete  was 
adequate  for  predicting  the  loading  response  of  a  concrete  structure. 

4.  The  model  developed  does  not  adequately  handle  unloading,  nor  does 
it  fit  in  detail  the  stress  state  dependence  of  yield.  The  latter  effect 
would  not  be  difficult  to  include  as  sufficient  data  become  available. 
Anisotropy  is  not  considered,  nor  are  environmental  effects  S"  as  temperature 
or  water,  or  preconditioning  effects  such  as  a  preshock  or  ..  preload. 
Time-dependent  effects,  creep  or  high-strain  rate  loading,  are  not  considered. 
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APPENDIX  A 


CONSTITUTIVE  EQUATION  SUBROUTINE 


SUBROUTINE  VARIABLE  DEFINITIONS 


AB1,  AB2,  AB3,  AB4,  AB4 
Al,  A2,  A3,  A4,  A5,  A6 
B1(I),  B2(I) 

AGO 

AKBULK 

AJ1I 

AJ2I 

C(I,J) 

DEHP 

DEPST(I) 

DF1 

DF2 

DS(I)‘ 

ELI 

EPSPI(I) 

FSI(I) 

FS2(J) 

FI 

F2 

HI 

H2 

PDE(J) 

SIGI(I) 

SIGDI(I) 

SI2P 

SI2PS 

SJ2 

TDS(I) 

TZL2 

Tl,  T2,  T3 ,  T4 ,  T5,  T6 

ZK1 

ZK2 

ZL1 

ZL2 


Terms  in  expression  for  Xj  and  x2 

Shear  modulus,  psi 

Bulk  modulus,  psi 

First  stress  invariant,  Oi,  psi 

Second  deviator  stress  invariant,  J2 

Matrix  of  elastic  constants 

1st  plastic  strain  invariant,  IjP 

Total  strain  increments 

Loading  increment  on  1st  yield  surface 

Loading  increment  on  2nd  yield  surface 

Stress  increment  associated  with  elastic  strain 

Lame  constant,  psi 

Plastic  strains 

Partial  derivative  vector  of  F]  with  respect  to  stress 
Partial  derivative  vector  of  F2  with  respect  to  stress 
1st  yield  surface 
2nd  yield  surface 

Strain  hardening  function  for  1st  yield  surface 

Strain  hardening  function  for  2nd  yield  surface 

Plastic  strain  increment 

Stresses 

Devi a  lor  stress 

2nd  plastic  deviator  strain  invariant,  /T2P 
Square  of  SI2P 
Square  root  of  o2,  JZ'2 

Stres-  increment  associated  with  plastic  strain 

Temporary  name  for  x2 

Terms  of  various  functions 

Constant  in  1st  yield  surface 

Constant  in  2nd  yield  surface 

Constant  in  associated  flow  law,  Xi 

Constant  in  associated  flow  law,  x2 


L 


% 


o  C)  o  o 


- THIS  SUBROUTINE  CALCULATES  stress  increments  from  an  input  of 

STPAIN  INCREMENTS,  STATE  OF  STRESS , PLASTIC  STRAIN,  AND  THE  YIELD 
FUNCTION  CONSTANTS 

- THIS  SUBROUTINE  USES  TWO  STRAIN*HARDENING  YIEIO  SURFACES 

SUOROUT INE  PLASTK(OEPST,SIGI,SI2P,EPSPI  ,ZK1,2K2> 

0 1  MENS  1  ON  SIGI  U)  .DEPST  <4>  ,0S<4>  ,SIGOI  <4>,C<4,<*>  ,FSi<4>  ,FS2(4)  , 

1  TDS  (  4  >  ,31(4)  , 32<4),EPSPI<4>  ,PDE<4) 

AGO=2 .  c.  +6 
AKBULK=3, £♦£> 

ELl=1.6666?E+6 
00  10  1=1,4 
00  1C  J=l,4 
10  C<I,J)=0. 

00  If  1  =  1,3 
00  15  J=1 , 3 
15  C ( I , J) =  EL  1 
DO  20  1=1,-+ 

23  C ( I , I) =C(I, I ) +2. *AGO 

AJ11=C!GI(1)+SIGI(2)+SIG!(3> 

A J2I  =  C  . 

uo  2'.;  i=i, j 

SIUUI(I)  =  SIGI  (D-AJll/3. 

2?  AJ21-AJ2I+SIG0I <I)**2 

AJ2I=,5*(AJ2I*Slf.I  (4)**2) 

SJ2  =  AJ2I**«  5 

C . CALCULATION  OF  ELASTIC  STRESS  INCREMENTS  FOLLOWS 

00  30  1=1,4 
3  0  0  S  < I >  =  0 . 

00  35  1=1,4 
00  35  J=1,h 

35  US(I)=0S(I)4C(I,J)*0£OST(J) 

ZL 1  =  0  # 

ZL  2  =  5  . 

T 1=  E  XP  < AJ1I/h030u. ) 

T2-1000.M12.2  -11.  ^Tl) 

T3=EXP((AJlI-  8 0 0 . ) / 7 0 0  . ) 
f  1  =  5 ■■J2-T2*  < 1 • -T  3)  -ZK1 
F?=.5*AJ2U2.*AJ1I*AJ1I-ZK2 

c - TEST  FCR  YIELDING  FOLLOWS 

IF (FI)  50,40,40 

,0  H1= (i.t-6) *  Cl.  ♦  100,*(“1.+EXP(120.*SI2P/<1.-,0001*AJ1I) )) )  / 

1  <?.-Z.E-5*AJlI) 

t*  =  . 275*3 1* < 1.-T3)  +  T2*T3/70C. 

DO  42  1=1,3 

42  FS1 (I)= (.5*5101 (I)-, l66bfc66?*AJlI)/5J2  +T4 
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rsi  u)  =  sic-i  (<*)  ♦.j/sjz 

5 J  IF (F2) 70,60,60 

60  T5  =  -<2.0E-14)*SJ2*(1.-EXP<-S.J2/(»Q00.>  >  *  £  XP  (- ( (  S  J2/945?.  )  *  *4 )  ) 
1  / (l.-(6.E“4>*AJlI> 

H2=-{7. 02 E-12)*(1,-EXPC(AJ1 1/23811.) *♦ 3 ) ) * E XP ( A J1I/60000 . ) +  T5 
DO  65  1=1,3 

65  FS2<T)=»5*SIGI(I)+3.6333333*AJ1I 
F  S  2 (4) =  . 5  *S IG I  (4) 

7 0  IF((F1.LT.  0.) . ANO. (F2.LT.  0.))  GO  TO  250 

C . CALCULATION  OF  PLASTIC  STRESS  INCREMENTS  FOLLOWS 

OO  02  1=1,4 

0 1  ( I)  =  0  • 

2  32  ( I )  =  0 • 

00  85  1=1,4 
no  85  i=i,4 

31(1)  01  ( I)  *C  (I  ,  J)  *FS1C  J) 

35  32 < I) =02 : I) +0 (I , J) *FS2( J) 

A  3 1  =  C  • 

A  R  2  =  0  , 

A  B  3  =  0  . 

A  3  4  -  G  « 

A  35=0. 

00  30  1=1,4 
A  3 1  =  A  3 1 +F  S l (1 >  *81(1) 

A32  =  A32-fFSl  (I)  *82(1) 

A33=A03+FS1 (I) *05(1! 

A3*=AB4+FS2<I)*82(I> 

9  J  A35  =  A65+FS 2 (I) *OS(I) 

A  1= 1.  +  2 ,*H1*A01 
A2  =  SJ2+2. *  H 1  *  A  3  2 
A3=2.*H1*A33 
Ah=. 3333333*H2» A32 
Ap=4.*AJlIf,3*33333*H2*AB4 
A  6  3333333*  H2*AUS 

I F <  F 1 )  91,92,02 

11  7L2  =  A6/A5 
GO  TO  <3€ 

9  2  IF  (1*2)  93,94,04 

93  ZL1=A3/Al 
GO  TO  9F 

^  T6= A1*A5-A2*A* 

7 L  1=  (A3*A5-  ■*. 2*. -JS)  / TC 
7l2=(A1*A6-m 3  *  A  h )  /  T  6 
?o  CONTINUE 
9  5  OF  1  =  0, 
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> 


OF  2  =  0 • 

00  luC  1=1,4 

TOS  <I>  =  DS (I) -81 <I) *ZL1-82<I>  *ZL2 
0F1=DF1+FS1 <I)*TOS(I> 

1)0  0F2=0F2+FS2(I)*T0S(I) 

T  ZL  2=  ZL  2 

IF { OF 1*ZL 1)  102,130,130 
102  IF ( 0F2*ZL 2)  105,110,110 
135  ZL  2  =  0 • 

GO  TO  95 
110  ZL 1=0 . 

ZL2=TZL2 
GO  TO  95 

130  IF<0F2*ZL2)  135,200,200 
135  ZL 2  =  0 • 

GO  TO  95 
230  CONTINUE 

C - INCREMENTING  OF  YIELO  SURFACES  FOLLOWS 

IF  (OF1*  ZL 1)  210,210,205 
235  Z*<l  =  ZKi+OFl 
210  IF<OF2*ZL2>  230,230,220 
220  ZK2=ZK2+DF2 
230  00  235  1  =  1,4 

SIGI(I)=SIGI(IMTDS(1> 

23  5  POEd)  =  FS1(I)¥ZL1+FS2<I)*ZL2 

0cHP=EPSPI(1)+£PSPI<2>+EPSPIC3) 

IF ( S I  2  P )  245,245,239 
239  00  240  1=1,3 

24  0  ni2P  =  SI2PM.5*EPSPim/SI2P-0£HP/6.>*PDEa> 

SI2P=SI2P+(,5*EPSPI(4)/SI2P) *POE  (4) 

GO  TO  248 

2-* 5  S 1 2PS=  (  (POC  (l)-POE  (2)  )  **2+ (P0E(2>-PDEC3n**2MPDE<3>-P0ECi>  >**2) 

1  /6.+P0E(4>**2 

S 1 2P=  SOPT (SI2PS) 

24d  DO  249  1=1,4 

2 49  E^SPI (I)=EPSPI (I)+P0E(I) 

GO  TO  3CC 

253  CONTINUE 

00  26 C  1=1,4 

250  SIGI(I>=SIGI\I)+OSCI) 

3)3  RETURN 

END 
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APPENDIX  B 


INELASTIC  THICK  WALLED  CYLINDER  PROGRAM 


Figure  04.  Inelastic  cylinder  program  flow  chart. 


£FC»,IS  MAIN, PAIN 

DIMENSION  ST  (AO  ,40)  ,  SR  (40,4Q)  ,SZ (40 ) ,ET (40 > , EP (40 )  »EZ(40)  • 

1  U(4C,4C) ,C(4, 4 ),F(40) ,X( 40) ,8(40) , A (40» 41)  ,PE (3) , 

2  S(4) , XP2 (40) #T (40,3) *  EE ( 3) ,A1(4) ,A2(4) ,GP1(3Q,3> ,GP2(30*3) , 

3  B  K 1  ( 3  0 ) »3K2(3Q)»D(3*3) 

33  FORMAT  (7(/,  11E12.4) ,/,El6.4) 

°E  AO  10.LTWC 
00  ICO  NT M= 1, LTMC 
pEAO  10,LSP,LTIM 
10  FORMAT  (215) 

H=2./LSP 

P=3.5 

PRINT  15,H,LTIM 

15  FORMAT  (//,5X,*THICK  WALLEO  CYLINDER*, 1 OX* /RADIAL  MESH  =  * , 
l  F8.4,*INCH*,10X,I5»*  LOADING  STEPS*)  * 

NN=LSP+1 

READ  25,S7(1) , E7IN, AP2,EIN 
25  FORMAT  (4F2C.8) 

PRINT  25, SZ (1) ,EZIN, AP2,EIN 

00  2C  1=1, NN 

ST(I,1)=0. 

SR(I,i)=0. 

ET(I) =EIM 
ER  ( I)  =  E IN 
EZ(I)=EZIM 
3<1 (I ) =C. 

0<2(I) =5.E+7 
VD2  (I )  =  A°2 

>1  U«I,1)= (R+H*(I-1))*EIM 
R(1)  =  C, 

point  3  G ,  (SP(I,1) ,I=1,NN) * (ST(I,1) ,I*1*NN) , <ERCI> »I  =  1*NN) , 

1  <5T (I) ,I  =  i,NN) , (XP2 ( I ) *  1  =  1 *NN) , (EZ(I) *  1= 1, NN) ,SZC1) 

ltim=ltim+i 

PEAD  40, (0(J) ,J=2,LTIM) 

°EAO  40  ,  (SZ(J) , J  =  2,LTIM) 

*0  FO°MAT  (8F10.5) 

00  37  J^2,LTIM 
SZ(J)=-S7(J) 

0(J)=  -P(J) 

37  SR ( 1, J) =B ( J) 

KK=3*NN-2 

S(4)=C, 

DO  IOC  J=  2 »  LT I H 
00  41  K=1,KK 
F ( <)=C» 
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00  41  L=1,KK 
'*1  A(K,L)  =  G. 

S(1)=SR(1, J-l) 

S ( ?) =ST  (1 ,  J- 1) 

^ ( 3) =SZ  C J-l > 

CP2=VP2(1) 

CK1=BK1 (1 ) 

CK2=BK2(1) 

CALL  PL  AST  ( 0 , S ,C<1 ,CK2 , CP ? , A  1 , A  2 , C) 

00  1  M=l,3 

GP1  (1  ,HM)  =A1  (MN) 

GP2 (1,KK) =A2(MM) 

1  T  <1,M)=C(  3,-M) 

CtttttFIRST  c  EQUATION 
P=0,B28*H 
A(ltl)=l. 

A(l,5)=-D (2,1) /P 
A (1,2) =0(2, l)/P-0(2,2)/R 

F(l)=ST  (1,  J-l)  -  (0(2,1)/P)*U(2,  J-l)4.(0(2,l)/P-0(2*2)/R)*U(l,  J-l) 

1  4-0(2,  ?>MSZ(J)-SZ(J-i)) 

C*****FIPST  b  equation 
a (2,5)30(1*1) /P 
A (2,2)= (0 (1,2)/R)-A(2,5) 

F  (2)  =  P(  J)  -  3  (J-l)  4-0  (1,1)  /P*  (U(2*J-1)-U(1*J-1))>  (0(1*2) /R)*U(1*  J-l) 
1  -0(1»3)MS7(J)-SZ(J-1)) 

CtttttFIRST  A  EON,  SECOND  POINT 
A(3»f)  =  l, 

A(3,3)=2.*H/(R4-H) 

F (3) =B( J) 

A (3,4) =-A (3, 3) 

NT=LSP-1 

C  ttt  yjfCENTPAL  A  EQUATIONS 
00  US  1=3, NT 
L  =  3*  (I-D+l 
M=L-1 

A(F,L+2)=1, 

A (N,l-4)=-l. 

A  (  M  ,  L  - 1)  =  2, *H/  (R4-(I-1)*H) 

U5  A(H,L)=-A(M,l-l) 

00  50  1=2, LSP 
S(1)=SR(I, J-l) 

S ( 2) =ST (I , J-l) 

S (3)=SZ (J-l) 

C°2=XF2(I) 

C<1=RK1(I) 
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OK2=PK2(I) 

cat.!.  PL  AS  T  <0,S,CK1,CK2,CP2,A1,A2,C> 

00  3  M M  =  i  ,  ?. 

0?1 (I =Ai{MM) 

0P2 (I ,MK> = A2(MK) 

3  T ( I , KM ) =C (3, MM) 

M=3*I-1 
L  =  M-1 

r**e**CENTRAL  0  EQUATIONS 
OT=9+  (1-1) *H 
A  CL,L-l)  =  i. 

A  (LfL  +  L)-=»3.5*0  ( 1 ,  1)  /H 
A<l»L-2)=-A<L»L4-M 
A (L,M)=-0(1,2)/PT 

F(L)=SR(I , J-l) -(3.5*0ti,l>/H>*CUCI*l,J-l)-U(I-l,J-l)>-0Cl,2>/RT 
1  *UCI,  J-l)  +0(1,3)  MSZC  J)-SZ(J-l)) 
cttt  **CENTPAL  C  EOUATIONS 
A  (M.,  L)  =  l. 

A  (M,l+4)=-3.5*0  (2,1)/H 
A  CM,  l-2)  =-A  (M,L  +  <*) 

A (M,M) =-0 (2,2)/PT 

E(M)=-0.5*O  (2,1)  /H*  (U(  I  +  l,  J-1)~UCI-1,  J-l) )  -0  (2,2>/RT*U(I,  J-l) 

1  +0<2,?) MSZCJJ-SZC  J-l>> 

50  continue 

L=3*LSP-3 
M  =  L  -  7 

C**«/*NEXT  TO  LAST  A  EQUATION,  NEXT  TO  LAST  POINT 
A (L»M)r-l. 

A  ( L  *  L  )  =  <R-MLSP-1)*H) 

A(L,L*i)=-A  (L,L) 

Zitttt NEXT  TO  LAST  3  EONCAOOITION  TO) 

A  (L  +  l  ,L+4)  =  -C'»5*0(i,l)/H 
A  <t*2,L+4>--<:.5*D<2f  1>/H 
L=5*LSP+1 
M  =  L  - 1 

PJ  =  R  +  ( L  Sr  *H ) 

C«f*<<LAST  A  EQUATION 
A  <M,M-3)  =  1. 

A  (M, V)  =  H/PT  *1.057 
Cttttt  tlkST  C  EQUATION 
S ( 1 )  =  SR  CNN  ,J-1) 

S ( 2 ) (UN  ,J-i) 

SC3) =SZ  ( J-l) 

CP2=XP2 (NN) 

C  K 1  =  F;  K 1  *  M  H ) 
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CK2=PK2 (NN) 

CALL  PLAST  ( 0 , S , CK1 ,CK2, CP2 , A 1 , 6 2, C) 

DO  4  MM=1,3 
G?1  (NN , MM ) = Al CMN) 

GP? (NN,  KM) =A2 (MM) 

4  T (NN  ,  MM) =C ( 3 »MM) 

P=i. G9l*H 
A  (L,M>=1. 

A (L  »L)  =  -  <0<2,2)/RT+0(2,l>  /P> 

A<L,L-2)=0(2,1)/P 

F(L)=-0  (2,  i)/o*<U(I,  J-1>-U(I-1,  J-i)  )  -0(2,2) /RT*U Cl,  J-1)*0C2,3)* 

1  (S7 CJ)-SZ CJ-1) > 

CALL  SOL ( NN , A , F , X) 

ST (  1 ,  J ) =X(1) 

U  Cl, J)=X(2> 

ST (NN , J ) -  X ( 3*NN- 3) 

U ( N  N , J ) =X (3*NN-2) 

00  6C  1=2, LSP 
L  =  7*I 

SP(I, J)=X CL-3) 

ST ( I , J) =X (L-2) 

S3  UCI,J) =X(L-i) 

SR (NM, J)=0. 

STCNN, J)=X(3*NN-3> 

UCNN,  J)  =  x(3*NN-2) 

00  7  5  1=1, MN 
FT(I)rU(I,J)/(R+<I-l)*H) 

00  71  1=2, LSP 

ER  (I)  =  (UlTfl, J) -UCI-1, J) )/C2.*H) 

7  1  EZ  (I)  =  EZ(  I)- (T<  1,1)*  CEP  CD- CUCI  +  1,  J-l) -UCI-i,  J-l)  >/C2.*H))*T  II,  2)* 
1  (ET (I) -U (T, J-l) /(R+ (I-1>*H)  )-CSZCJ)-S7CJ-l>>>  /TCI,3> 

ER  (1)  =  (UC2,  J>  -UC1,  J)  >/H 

EZ(1)=£Z(1> -CT(1,1)* CER(1)-CUC2, J-l) -U Cl, J-l) > /H) ♦! Cl, 2) * CET C i> - 
1  U  ( 1 , J-l ) /R) -  CSZ(J)-SZCJ-1)))/TC1,3) 

ER (NM) = (U(NN, J) -U CNN-1, J) ) /H 

EZ (NM) =  EZ (NN) -(T (NN,1>*  CER CNN)- CU CNN, J-l>-U CNN-1, J-l) )/H)  +  TCWN,2) 

1  * (ET (NN)-U(NN, J-l) /(R+LSP*H))-CSZC J)-SZC J- 1) ) ) /T CNN, 3) 

S ( 3) =SZ  (J) 

CALL  PROP  CC ) 

00  90  1=1, NN 
00  31  M  =  1 , 3 
31  EE  CM) =e. 

S(i)=SR(I,J> 

S ( 2 )  =  ST  (I, J) 

00  82  M=1 , 3 
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00  8  2  L-lfS 

82  EE<K)=£E<M>+C(M,L)*S<L> 

P£ (1)=EP<I>-EEC1> 

P£ <2>=ET<I> -EEC2) 

?E(3)=£Z(I)-££(3> 

PM  =  (PE  (1)  +PE(2>  «-PE  (3  J>/3. 
no  83  K  =  l,3 

83  P"  (Ml  =PE(M)  -PM 
X?2(I)=0. 

00  8i*  M=1  *  3 

84  XP?(I)=XP2(I»+PE(M) *PE(H) 

93  YP2(I) =  SQRT (XP2 (I) *0.5) 

00  22C  1=1, NM 

OBKl=GPltI,l>MS°<I»J>-SRCI,  J-l) >*GP1CI,2> MST  CI» J)-ST Cl,  J-l)) 
1  GP1(I,3)MSZ(J)-SZCJ-1>) 

00<2=G?2(I, 1) * (SRC I, J)-SPC I , J-i) ) ♦GP2  C 1 , 2) * ( ST  Cl , J) -ST Cl , J-l) > 
1  GP2(I,3)*(SZCJ)-SZCJ-1>> 

TP ( C°K1)  210*210,205 
235  BK1  (T)-P<ltlU08Kl 
210  IF  <  Qn<2)  220,220,215 
215  OK2(T)=GK2tI)+OPK2 
0  CONTINUE 

PRINT  3^.,  A=1,NM)  ,CSTCI,J)  ,I=1,NN),  CERCI)  ,Isl,NM), 

1  <ET  ilf  r  ,  1  =  1 , NN)  ,  CE Z  ( I )  *1-1* NN)  ,  (U  Cl,  J)  ,  1=1  ,N) 

2  S7(J> 

130  CONTINUE 

PNO 


(for, is  supj, vu‘u 

S:J3°CUTIMc  po0P(C) 
DIHEMSIOM  C<4,4) 
EM- 4. 909091E+E 
EU=0. 2272727 
no  5  1=1,3 
00  5  J=  1, 3 
5  C (I, J)=-EU/EM 

no  ic  l=i,3 

10  C(I,T)=1./EM 
Pc  TUPN 
ENO 


CEO*,  IS  SIJP2.SUB2 

SUBROUTINE  SOL ( NN , A , F , X ) 

DIMENSION  A  (40,41)  ,F  (40) , X ( 4 0 ) 

N=.3*NN-2 

NP=N+1 

OO  2  1=1, N 

2  MI,NP)=F(I> 

00  9  I=1,M 
IP-1+1 

00  9  J=1,N 
IF  <I-J>4,9,4 
4  G  =  -A ( Jt I) /ft  (1 ,1 > 
no  3  K  = I P , N  P 

3  A(J,K)=A< J, K) +G*A (I, K) 

9  CONTINUE 

no  ic  i=i, m 

13  XCI)=A(I,NP)/ACI,I) 

°E  TUPM 
END 
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(FOR, IS  SU°3,SU03 

SUBROUTINE  PLAST(0,SIGI,Z<l,ZK2,?I2P,FSi,FS2,C> 

C - THIS  SUBROUTINE  CALCULATES  STRESS  INCREMENTS  FROM  AN  INPUT  OF 

C  STRAIN  INCREMENTS,  STATE  OF  STRESS , PLASTIC  STRAIN,  AND  THE  YIELD 

C  FUNCTION  CONSTANTS 

0 - THIS  SUBROUTINE  USES  TWO  STRAIN*H6R0ENING  YIELO  SURFACES 

DIMENSION  SIGIU),  SIGOI  (4)  ,£(4,4)  ,FS1  (4)  ,  FS2<4>  , 

1  01(4), 32(4),  Pl(3,3> ,P2(3,3), 

2  P3(3,3) , PA (3,3) ,OC3,3) 

AG0=2 , E+6 

A  K3ULK=  3. E  +  fc 
ELi=i.66667E*6 
00  1C  1=1,4 
FSl <I)=C. 

FS2(I)=C. 

00  1C  J=1 » 4 
10  C(I,J)  =  C. 

00  IE  1=1,3 
00  IE  J=i »  3 
1?  C  ( I ,  J)  =EL  1 
00  2C  1=1,4 

25  C(I,I)=C(I, I)+2.*AG0 

AJH=SIGI  (D+SIGK2)  +SIGI  (3) 

A  J  2 1  =  0  • 

DO  2E  1=1,3 

rIF,ni  ( I )  =SIGI  (D-AJ1I/3, 

AJ2I=AJ2T+SIG0I (I)**2 
AJ2I=.5*(AJ2I+SIGI(4>**2) 

SJ2=A J2I**,5 
7L 1 =5 • 

ZL  ?  =  C  • 

Tl==y?< AJ1I/40000.) 

T  2= 1C  G  G  »  * (12,2  -ll.*Ti) 

T3=E7P( (AJ1I-  303, )  /70G • ) 

F1=SJ2-T2*(1.-T3)-ZK1 

F2=»E*AJ2I+2,*AJ1I*AJ1I-ZK2 

r - TEST  FOP  YIELDING  FOLLOWS 

IF(F1)  E0,4C,4Q 

4Q  hi=  (i.E-6) * (1,  ♦  10O.*(-l.*EXP{i23.*SI2P7(l.-.OOOi*AJlI> >)>  / 

1  ( 2, -7 ,E“5¥ A Jll ) 

T4=  ,27S*TlMi.-T3)  *  T2*T3/700. 

00  u2  1=1,3 

4?  FSl (T)= (,5*SIGI (I)-. 16666667*AJ1I) /SJ2  +T4 
FS1(4)=SIGI(4)*,5/SJ2 
50  IF  (F2)7C,60,e0 

iG  T5=  -(2,0E-14)*SJ2*(l.-EXP(-SJ2/400O.))  *  £XP (- ( (S J2/9457. > **4> ) 
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1  /(1.-(F.E-4)*AJ1I> 

H?=- I7.02E-12)* <1.-EXP<  < A J1 1/2 3811.)** 3) >*EXP<A J1 1/6 0000.) *T5 
no  65  1=1,3 

85  FS2(I>=„5*SIGI<I>*3.8333333*AJ1I 
FS2<4)=,S*SIGI(4) 

r:  ifkfi.lt.  o.  >. and. (F2.lt.  o.n  go  to  250 

- CALCULATION  OF  PLASTIC  STRESS  INCREMENTS  FOLLOWS 

DO  82  I  =  i,4 
31 ( I )  =  0 » 

82  fl2(I) *C. 

00  85  1=1,4 
00  86  J=l,4 

3i(i)=eim  +cd,  j)*fsk  j) 

35  e?(I)=B2(n+C(I,  JMFS2CJ) 

AB1=C . 

A92=C . 

A93=C . 

A04=r. 

A B5  =  0  • 
r'0  9 C  1  =  1,4 
A91=C31+FS1(I)*31(I) 

A9?=AB2*FSl (I)*R2(I) 

AB*,  =  AB4+FS2(I)*B2(I) 

93  CONTINUE 

Al=l.+2.*H1*A91 
A2=SJ?+2. *H l*  A92 
A4=. 3333333*H2*A82 
A  5  =  4. *AJ1I  +  .3  333333*H2*A84 
00  28  C  1  =  1,3 
00  2fiC  J=  1 ,  3 

pki,  j)=9i(i)*ni(j> 

02(1, J)=3l(I)*82(J) 

P4(I, J)=B2(I) *91<J) 

233  *>3(I,  J)=R2(I>*32(JJ 
A  9=  2 ,  * H 1 
A 9=  0 «  33 33 3* H 2 
29Q  A io=A l*A5-A2#A4 
00  2«5  1=1,7 
00  295  J=  1 ,  :< 

295  C(I,J)=  C <1 , J)- (PI (I , J) *A5* A8-P2 (I, J> * ( A2* A9)  *P3 < I , J) *A1*A 9 

1  -P4  (I  ,  J> *A4*A8) /A10 

25  3  CONTINUE 
338  CONTI NUF 

00  310  1=1,2 
0(I,3)=C(I,3)/Ci3,3) 

00  31  0  J=  1 , 2 

310  DU,  J)  =C(I,  J)  -Cl  I,  3)  *C(  J,3)  /C(3,3> 

PE  TUPN 
ENO 


1  .6 


nr  nr.  n  wai  nrrao>:  *-<wt  '<  r.T:  ’)•>' 


APPENDIX  C 

DERIVATION  OF  EQUATIONS  FOR  A]  AND  a2 

In  Section  V  it  was  stated  that  equations  13  and  15  could  be  developed 
by  substitution  of  previous  equations.  This  will  be  shown  in  this  Appendix. 
Starting  with  the  hardening  rule  assumed  in  the  form  of  equation  10, 


Section  V 


d^P  =  hj  dKj 


where  by  d/T'p  it  is  meant 


d/Hp  .  j  (dc,p-dc2P)2  +  (^2p-dt3p)2  +  (dc3p-dtjp)2 


(dL  12)2  +  (dfzP);'  +  (d£31)2 


As  explained  in  Section  V,  since  d/[)p  involves  strains  from  both  yield 
surfaces,  the  flow  rule  given  by 


d'  iiP  "  ;jFi  Ai  +  A2 


should  be  considered.  This  equation  (7)  is  substituted  into  (10)  to  obtain 
the  desired  result.  Making  this  substitution 


■jF  ]  >  1  +  •; 


■JF ;  A]  +  jF2  >2\] 


+  r  *  f  .  > .  +  i  f  \ .  *) + 

;  i  / 

-  ’  ’  i :  -  *'■  i  r 


!=  h.  dK, 


Since  the  yield  surface  segments  are  functions  of  the  stress  invariants 
/T)  and  only,  the  derivatives  can  be  written  as 


3oij 


aFj  s/Jj  +  aFj  aO i 
3*^2  3o ^ j  3 J i  3o^j 


(48) 


and  similarly  for  F2.  Noting  the  definition  of  the  stress  invariants,  the 
derivatives  can  be  written  as 


bo..  2/J2  6/J2 

ij 


and 


dJ  1 


6ij 


dO  .  . 

IJ 


where 


ij 


0  if  i  f  j 
1  if  i  =  j 


(49) 


substituting  equation  C-3  into  C-2  and  rearranging  gives 


aFi  j 

- 

3^2 

V1 

+  dF2 

•  3  >'  J  2 

Vaon 

3a22 

) 

d/$2 

3j  j 

-  3  J  ]  \ 

1  X1  + 

cj  F  2  j 

pj, 

;j0  1  1 

d022  j 

3  j2 

\^U 

d'/Jn  \  ^(7  ]  ]  9^22 

0|_V: 


+  T/3F,  3/J)  +  DF-.oJj  ^  X  j  +  I  0^2  3  +  3F2  ^Ar  1  2+ 


l\^2  j  9oi2 

=  hi  dK, 


( 3  /J'2  do  j  2  9  J  J  do  • 


and  using  49  simplifies  this  to 
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j  1  ~  e>  F  i  /  Q  1  1  ?  2  \  A1  +  aF2  /an~a22  \  A2~]  2  +  " 

/  6  Lj/J2  V  2^2  /  a  ^  \  2^  /  _ 


+  roFt  a  12  Aj  +  3p2  OJ2  a2*|^  +  “"’( 

2/T2r  a/J^  2/^2  J  ’ 


=  h,  d K j 


Factoring  the  terms  gives 


2JJ'2  )  (  oQ 


j F j  Xj  +  3F2  X2  )  1  (cjj-Ojj)2  + 


-  +  o122  + 


hi  dKj 


Noting  that  the  last  term  is  just  /J’ 


1  j  jF,  >.j  +  ;.-F2  x2  )  /J2  =  hj  dKj 


2/J-  / 


and  finally  equation  13  of  Section  V  is  obtained  as 


3Fi  >,  +  OF,  =  2  h,  dKj 


Equation  15  is  obtained  in  a  similar  manner.  Starting  with  the 
hardening  rule  assumed  in  the  form 


dr  P  =  h2  dK2 
VC  ^ 


where  bydyrp  it  is  meant 


d>  vcP  =  (d  np  +  d •. 2 2 p  +  dr35p)  segment  F? 


As  discussed  in  Section  V,  equation  31  involves  the  compaction  strain 
from  the  F2  segment  of  the  yield  surface  only.  Thus  the  strains  obtained 
from  this  yield  segment  only  should  be  used.  These  are  given  by 


(deij  )p  =  ^  A2 

2  3c.. 

1J 

Substituting  equation  14  into  equation  11  gives 


(14) 


3F  - 


+  dF p  X ,  +  3  F  ->  X  =  h7  d  K  - 


(53) 


aC  i  ] 


oc22 


°  0  3  3 


Since  F2  has  been  assumed  to  be  a  function  of  and  di  only,  the 
derivatives  can  be  expressed  as  done  previously  in  C-3  for  F}.  Making 
this  substitution  gives 


3F2 

cU+c22+a33  '  3J] 

+  3F2 

[3]1 

-  /T* 

j  v  J  2 

L  2/J'  6/0),  J 

^7 

=  h2  d«2 


(54) 


Noting  that  the  first  term  is  zero  gives 


3  3F2  X2  -  h2  dK2 
3  J , 


(55) 


and  finally  equation  15  of  Section  V  is  obtained  as 


3F„  x,  =  1/3  h,  dK„ 

/  2  2  c 


(15) 
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